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An alternative method is suggested for the description of the velocity and pressure fields 
in an unbounded incompressible viscous fluid induced by an arbitrary number of spheres 
moving and rotating in it. Within the framework of this approach, we obtain the general 
relations for forces and torques exerted by the fluid on the spheres. The behavior of the 
translational, rotational, and coupled friction and mobility tensors in various frequency 
domains are analyzed up to the terms of the third order in the dimensionless parameter 
b equal to the ratio of a typical radius of a sphere to the penetration depth of transverse 
waves and a certain power of the dimensionless parameter a equal to the ratio of a typical 
radius of a sphere to the distance between the centers of two spheres. We establish that 
the retardation effects can essentially affect the character of the hydrodynamic interactions 
between the spheres. 
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I. INTRODUCTION 



The theory of hydrodynamic interactions between particles of finite size immersed in a 
viscous fluid is of great interest for numerous problems of the modern physics and tech- 
nology such as the theory of Brownian motion [1-6] stimulated by the results of numerical 
calculations on computers (including "Brownian dynamics" in computational physics [7]), 
sedimentation [8,9], etc. As a rule, the central problem of the theory is to determine forces 
and torques exerted by the fluid on particles moving and rotating in it with given trans- 
lational and angular velocities. To this end, the linearized Navier-Stokes equation for the 
fluid is usually used [8,10]. Starting from this equation, a great number of rigorous results 
has been obtained for a single sphere in a fluid [11-16]. For more than one sphere, there 
are only a few rigorous results obtained for two spheres in some particular cases of their 
stationary motion and rotation [17-20]. As for the nonstationary case, any exact results for 
the fluid velocity and pressure induced by two and more interacting spheres are absent. 

The situation with finding approximate solutions of the linearized Navier-Stokes equation 
is more advanced (first of all, at the expense of the classical method of reflections proposed 
by Smoluchowski and its further modifications). The comprehensive review of the results 
obtained in this way is given in [8] including the cases of mixed slip-stick boundary conditions 
at the surfaces of the spheres [21-25] and permeable spheres [26-28]. It follows that the 
majority of these results have been obtained in the stationary case and mainly for two 
spheres, and the possibility of an application of the corresponding methods to nonstationary 
many-sphere problems is not clear. 

In [4,29], on the basis of the method of induced forces proposed in [30], it was developed 
a procedure of the derivation of a hierarchy of equations for irreducible force multipoles 
induced in spheres immersed in a fluid. Using this procedure, the friction and mobility 
tensors of a system of an arbitrary number of spheres are obtained in the stationary case 
and the results are represented in the form of power series expansions in the dimensionless 
parameter a equal to the ratio of a typical radius of a sphere to the distance between the 
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centers of two spheres retaining in the series terms up to a certain order in this parameter. 
This method essentially differs from the methods of reflections and gives one the possibility 
to find the friction and mobility tensors without knowledge of the fluid velocity and pressure 
in an explicit form. Furthermore, the method is free from the restriction to the stationarity 
of the system. Its generalization to the nonstationary case has been done in [31] where 
the expressions for the mobility tensors have been represented in the form of power series 
expansions in two dimensionless parameters: a, which is typical of the stationary case, and 
the parameter b, which is the ratio of a typical radius of a sphere to the penetration depth of 
transverse waves. At finite frequencies, the expressions for the mobility tensors are obtained 
in [31] up to the third order in these two dimensionless parameters, namely, up to terms of 
the order of b n a m with n + m = 3, where n, m — 0, 1, 2, 3. 

There are some other methods [32-38] proposed for determination of an approximate 
solution of problems of many-sphere hydrodynamic interactions in a fluid in both the sta- 
tionary and nonstationary cases. As compared with the method of induced forces used in 
[4,29,31], these methods enable one to calculate not only forces and torques exerted by the 
fluid on spheres but also the velocity field of the fluid induced by these spheres (since the 
latter problem is much more complicated, usually only the hydrodynamic forces and torques 
are calculated.) 

The original analytic-numerical method has been used in [36] to study the quasistatic 
hydrodynamic interactions in a system of N spherical particles. Unlike the other methods, 
this method is not based on using the conventional scheme of the representation of the 
quantities under study as series expansions in powers of the parameter a but reduces the 
problem to the specific solution of an infinite system of linear equations. With this method, 
in [36], the stationary mobility tensors for a system of two interacting spheres of equal 
radii have been obtained and analyzed in detail. In [37], the method has been extended to 
the description of the retarded hydrodynamic interaction between two equal spheres in an 
unbounded fluid. The nonstationary mobility tensors of this system have been calculated 
and plotted versus the variable b for several values of the parameter a. Note that these 
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numerical results are not restricted to the range of small values of b like in [31] where the 
corresponding analytic relations for the mobility tensors are deduced. 

As for the method developed in [34,35], it should be noted that the final results presented 
in these papers are expressed in terms of multiindex tensors defined as tensors inverse to 
certain rather complicated multiindex tensors but the explicit expressions for these inverse 
tensors are not given. Naturally, without knowledge of the explicit form for the inverse 
tensors, the expressions given in [34,35] cannot be reduced to the known results obtained 
by other methods. Moreover, it can be shown [39] that one of the inverse tensors defining 
in [34,35] the relations for the forces and torques exerted by the fluid on the spheres and 
the velocity field of the fluid induced by the spheres does not exist (even for noninteracting 
spheres) because the corresponding direct tensor is singular. In [39], the reasons for this 
fact have been analyzed in detail, and the procedure has been proposed for finding the time- 
dependent distributions of the velocity and pressure fields in an unbounded incompressible 
viscous fluid in an arbitrary external force field induced by any number of spheres immersed 
in the fluid as well the forces and torques exerted by the fluid on the spheres. In the 
stationary case, this program has been performed for the derivation of the general relations 
for the velocity and pressure fields of the fluid, the friction and mobility tensors for any 
number of spheres in a fluid in given external force fields and these results are represented in 
the form of power series expansions in the dimensionless parameter a retaining in the series 
terms of the second iteration inclusive. 

However, it is well known [14,16] that even in the case of a single sphere an account of 
time dependence can lead to significant changes both in the fluid velocity induced by the 
sphere moving in the fluid with time-dependent velocity and in the force exerted by the fluid 
on the sphere, namely, the asymptotics of the fluid velocity, being of the order of 1/r, where 
r is the distance from the center of the sphere to the point of observation, for the stationary 
motion of the sphere, is replaced by 1/r 3 and the time-dependent force exerted by the fluid 
on the sphere, in addition to the ordinary Stokes friction term, contains a term depending on 
the acceleration of the sphere and a memory term associated with hydrodynamic retardation 
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effects. In the nonstationary case, hydrodynamic interactions between spheres also change, 
which is especially significant for large distances between the spheres. For example, certain 
mobility tensors, which are proportional to the inverse distance between spheres to certain 
powers in the stationary case, exponentially vanish with increase in the distance between the 
spheres in the nonstationary case [31]. The account of the main frequency-dependent terms 
in the mobilities in the case of small frequencies leads to long-time tails for the velocity 
correlation functions of two different spheres [31]. 

Naturally, it is of interest to investigate the nonstationary velocity and pressure fields of 
the fluid induced by several spheres immersed in the fluid. However, this problem cannot 
be solved with the use of the methods exploited in the pioneering papers [4,29,31]. It is also 
of interest to extend the analytic results for the nonstationary mobilities obtained in [31] 
retaining higher-order terms as well as to derive frequency-dependent friction tensors for a 
system of spheres. These problems are considered in the present paper. 

In Sec. 2, we give the general relations for the fluid velocity and pressure induced by 
a system of arbitrary number of spheres and the total forces and torques exerted by the 
fluid and external force fields on the spheres. The results are obtained in the explicit form 
without imposing any additional restrictions on the size of spheres, distances between them, 
and frequency range. These relations are expressed in terms of the unknown induced surface 
forces distributed over the surfaces of the spheres. For the determination of the unknown 
harmonics of the induced surface force densities, a closed infinite system of linear algebraic 
equations is given. On the basis of analysis of the quantities contained in this system, we 
propose to seek its solution by the method of successive approximations using the system 
of noninteracting spheres as a zero approximation and formulate a rule for determination of 
the fluid velocity and pressure fields as well as the forces and torques exerted by the fluid 
on spheres up to terms of b 3 and a given power of the parameter a. 

In Sec. 3, in the approximation of noninteracting spheres, we obtain the fluid velocity 
and pressure induced by a system of spheres and analyze their behavior in the far zone. The 
relations for harmonics of induced surface force densities obtained for the first and second 
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iterations are given in Sees. 4 and 5, respectively. We derive the general relations for the 
forces and torques exerted by the fluid on the spheres up to the terms of the order of b 3 and 
certain powers of the parameter a corresponding to the second iteration. We analyze the 
friction tensors in the particular cases: for sufficiently large distances between the spheres 
and in a low-frequency range. We show that the account of retardation effects essentially 
changes the character of hydrodynamic interactions between the spheres. 

In Sec. 6, the general relations for the translational, rotational, and coupled mobilities 
are derived up to the terms of the order b 3 and a 2 , a 4 , and a 3 , respectively. We estimate 
contributions of three-sphere interactions to the mobility tensors and show the possibility of 
coupling between translational and rotational motions of the same sphere. We investigate as 
well the influence of nonstationarity on the behavior of the mobilities in the low-frequency 
range and in the case of large distances between the spheres. 



We consider N homogeneous macroscopic spheres of radii a a and masses m a , where 
a = 1, 2, . . . , N, immersed in an unbounded incompressible viscous fluid at rest at infinity 
in an external force field and determine the fluid velocity and pressure induced by these 
spheres as well as the forces and torques exerted by the fluid on the spheres. 

In the linear approximation, the fluid is described by the linearized Navier-Stokes equa- 
tion 



II. GENERAL RELATIONS 




(2.1) 



and the continuity equation 




(2.2) 



Here, P{r, t) is the stress tensor of the fluid with the components 




(2.3) 



6 



p(r,t) and v(r,t) are the pressure and velocity fields of the fluid, p and rj are its density 
and viscosity, respectively, 5^ is the Kronecker symbol, and F ext (r,t) is the external force 
acting on a unit volume of the fluid. In the present paper, for simplicity, we consider the 
case of a conservative external force field, i.e., 

F ext (r,t) = -pV<p(r,t), (2.4) 

where ip(r,t) is the potential of this force field. (The corresponding results for the general 
case where the external force field has both potential and solenoidal components are given 
in [39]) 

The position of the center of sphere a at the time t is defined by the radius vector 
R a (t) relative to the fixed Cartesian coordinate system with origin at the point O (in what 
follows, the system O). In parallel with this fixed system, we also introduce N local moving 
Cartesian coordinate systems with origins O a at the centers of spheres coinciding with their 
centers of mass (in what follows, systems O a ) so that sphere a does not move relative to 
system O a . Therefore, the radius vector r of any point of the space can be represented in 
the form r = R a + r Q , where r a is the radius vector of this point relative to the local system 
O a . The motion of the spheres is described by the equations 

™«^r = Fm (2 - 5) 

I*^-=T?(t), (2.6) 

where 

U a (t) = (2.7) 

and f2 a (i) are, respectively, the translational and angular velocities of sphere a, I a = 
(2/h)m a a 2 a is its moment of inertia, F*"*^) is the total force and T^°'(t) is the total torque 
(here and below, all torques corresponding to sphere a are considered relative to its center) 
acting on sphere a 
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F t : t (t) = F e :\t) + Fl(t), (2.8) 
T?(t)=TZ t {t)+T' a (t), (2.9) 

where F ext {t) and T ext {t) are the force and torque acting on sphere a due to the external 
force field and F* a {t) and T* a {t) are the force and torque exerted by the fluid on sphere a 



Fi(t) = - J P(r,t)-n a dS a , (2.10) 

S a (t) 

T{(t) = - j ((r-R a (t))xP(r,t))-n a dS a , (2.11) 

S a (t) 

where S a (t) is the surface of sphere a at the time t and n a is the outward unit vector to 
this surface. 

Equations (2.1)-(2.3) describing the motion of the fluid are defined in the domain r a > 
a a , a — 1, 2 . . . , N, where r a = \r a \ = \r — R a (t)\. Following [39], we extend the external 
force field F ext (r,t) given in Eq. (2.1) for r a > a a to the domains r a < a a so that the 
extended quantity is described by the same analytic expression in the entire space as in the 
domain r a > a a (the analytic extension) and introduce new functions v(r,t) and p(r,t) 
that coincide with the fluid velocity and pressure in the domain r a > a a being defined in 
the entire space as follows: 



N 



a=l 



1 - Q(r a - a a ) 



U a (r,t)\, (2.12) 



p(r,t) = f;|e(r a - a a )p(r,t) + [l - Q(r a - a a )]{pW(r,t) - p (r ■ ^^)}}, (2-13) 
where 

p(v)(r,t) = - PV (r,t)+p mf (2.14) 

is the pressure of the unbounded fluid in the external force field F ext (r,t) in the absence of 
spheres and the constant p m ^ is the fluid pressure at infinity where the external force field 
is absent, Q(x) is the Heaviside function, and 

U a (r,t) = U a (t) + (fi a (f) x r a ). (2.15) 



As a result, the quantities v(r,t) and p(r,t) satisfy the continuity equation (2.2) and 
the equation 

p dv ^ + div P(r, t) = F ext (r, t) + F md (r, t) (2.16) 

given in the entire space, where the quantity Pij(r,t) coinciding with the stress tensor for 
r a > a a is given in the entire space by relation (2.3) with v{r, t) and p(r, t) instead of v(r, t) 
and p(r,t), respectively; the quantity F ext (r,t) is given in the entire space by the analytic 
expression for the external force field, and the induced force density F md (r,t) has the form 

N 

F ind (r,t) = J2 F a d (r,t), (2.17) 

a=l 

where F™ d (r, t) is the induced force density for sphere a containing the volume, F^ md (r, t), 
and surface, F^ md (r,t), components 

F™ d (r, t) = F<y) ind {r, t) + F ( a s)ind (r, t) (2.18) 

given inside the volume V a = (4/3) na? a occupied by the sphere a and on its surface S a (t), 
respectively, and defined as follows [39]: 

F<P ind (r,t) = &(a a -r a )p^- t (n a (t) x r a ), (2.19) 
F^ ind (r,t) = jdn a 5(r - R a (t) - a a )f a (a a ,t), (2.20) 

where f a (a a ,t) is a certain unknown density of the induced surface force distributed over 
the surface S a (t) and a a = (a a , a , <p a ) is the vector directed from the center of sphere a 
to the point on its surface characterized by the polar 9 a and azimuth ip a angles in the local 
spherical coordinate system O a , dfl a — sin 9 a d9 a dip a is the solid angle. 

Force (2.10) and torque (2.11) exerted by the fluid on sphere a are represented as follows: 

Fl(t) = F a (t) + FZ(t)-P?(t), (2.21) 
Ti(t)=T a (t)-fT(t), (2.22) 



where F a (t) and T a (t) are, respectively, the force and the torque exerted on sphere a due 
to the induced force density distributed over its surface 
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F a (t) = - JdrF^ md {r,t) = - JdQ a f a (a a ,t), (2.23) 

T a (t) = -fdr (r Q x F^ md (r,t)) = - jdil a (o Q x f a (a a ,t)), (2.24) 

F^it) and T^it) are, respectively, the force and the torque acting on a fluid sphere 
occupying the volume V a instead of sphere a due the external force field analytically extended 
to this domain, 

F™\t) = JdrF ext (r,t), (2.25) 

??(t) = Jdr (r a xF ext (r,t)), (2.26) 



and 

K® = (2.27) 

is the inertial force that is necessary to be applied to a fluid sphere of volume V a to ensure 
its movement with the acceleration dU a (t)/dt [30], where m a = pV a is the mass of the fluid 
displaced by sphere a. 

In what follows, we expand various quantities in the Fourier integral with respect to the 
time of the form 

1 °° 

A(r, t) = — ( duo A(r, u) exp(-iut). (2.28) 

27T J 

— oo 

Within the framework of the linear approximation both with respect to the fluid velocity 
and the velocities of the spheres [31], the solution of problem (2.2), (2.3), (2.16) for the 
required Fourier transforms of the quantities v(r,u) and p(r,u) can be represented as 

N 

v(r,u) = Y,vt d (r^), (2-29) 

(3=1 

N 

p(r,u)=p (0 \r,u)+p( md \r,ij)=p( \r,uj) + Y / P7 d (r,L;), (2.30) 

0=1 

where 

<(r, U ) = v^ md (r, u) + vf md (r, w ), (2.31) 
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{V)ind ( x ., ,2 2 1 1 A1 exp(ik- (r - R p )) j 2 {Ka p ) , x 



! 3?7 (2tt) 3 J k 2 + 

—oo 

1 00 

v%'"~{r,u) = j—^ J dk exp (ik • (r - Rp)) JdVtp exp (-ifc • ap) S(k,uj) ■ fp(ap,u) 



(S)ind, 



(2.33) 

pf(r,a;) = -^3s J rffc ~ ~ 'p ~ ~ fop exp(-ik> a p )(k> f^u)). (2.34) 

— oo 

Here, k = ^Ju]Jiu) = (1 — i sign a>) /<5, 5 = ^2v/\uj\, v = rj/p is the kinematic viscosity of 
the fluid, £0 = Qitrjdp is the Stokes friction coefficient for a sphere of radius ag uniformly 
moving in the fluid, &g = /tag is the dimensionless parameter that characterizes the ratio 
of the radius dp of sphere f3 to the depth of penetration 5 of a plane transverse wave of 
frequency to created by an oscillating solid surface into the fluid [14], (\bp\ = \/2ap/5), j n (x) 
is the spherical Bessel function of order n, and 

S(k,u) = ^^(1 ~ n k n k ) (2.35) 

is the dynamic Oseen tensor (J is the unit tensor and n k = k/k is the unit vector directed 
along the vector k). 

Note that Eqs. (2.29)-(2.34) are valid for any r. For r a > a a , the quantities v^ md (r, u) 
and p l p d (r,u)) can be interpreted, respectively, as the fluid velocity and pressure at the 
point r generated by the induced surface force F^ md (r,uj) distributed over the surface of 
sphere (3, whereas the quantity Vp V ^ nd (r, 00) has the sense of the fluid velocity at the point 
r generated by the induced volume force F^ md (r,u) given in the volume Vp occupied by 
sphere (3. 

It is convenient to represent the quantities Vp V ^ nd (r,uj), v^ md (r, u), and p l p d (r,uj) as 
the following expansions in the spherical harmonics: 

v ( p md (R a + r a , u) = 2a/tF v %im d ( R cn r ^ uj)Yi m (9 a , <p a ), (2.36) 

lm 

vf )md (R a + r a , u) = 2v^r ^ Vp S J™ d (R a , r a , Lv)Y lm (6 a , ip a ), (2.37) 

lm 

p i p nd (R a + r a ,u) = 2V*^P$ m (R a ,r a ,u)Y lm (6 a ,<p a ). (2.38) 



lm 
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Here, the spherical harmonics Yi rn (9,(p) are defined as [40] 

Y lm (6, <p) = 6XP ^ V?) e Zm (cos^), I = 0, 1, . . . , -Z < m < Z, 

V Z7T 



m—\m\ 



2/ + 1 (/ - Iml)! i ?77 1 / \ , 

^— WrT' w ' (2 - 39) 



oo I 

and P{ n (x) is the associated Legendre polynomial. For short, the symbol J2 means J2 ■ 

lm 1=0 m=— I 

As was mentioned above, the position of any point of the space defined by the radius vector 
r relative to the fixed system O can be represented in the form r = R a + r a , where R a is 
the radius vector that defines the position of the center of certain sphere a and r a is the 
radius vector directed from this center to the point at hand. We take the quantity r a such 
that it corresponds to the minimum difference — ag, where (5 — 1, 2, . . . , N, if the point of 
observation belongs to the domain occupied by the fluid, which means the minimum distance 
from the point at hand to the surface of sphere a. If the point of observation lies inside 
the domain occupied by sphere a, then r a is the radius vector directed from the center of 
this sphere to this point; r a = (r a , 9 a , ip a ), where 6 a and ip a are, respectively, the polar and 
azimuth angles of the vector r a in the local spherical coordinate system O a . 
The expansion coefficients in the series (2.36)-(2.38) have the form 

x^ m2 (e a/3 ,$ a/3 ), (2.40) 
vfjZ^Ra^u;) = £ T l ^l(r a ,a,,R a ,,u) ■ f^u), (2.41) 

pgf imi (R a , r a , w ) = £ D\\Z\ (r«, a p , R a p) • f p , hm2 (w) (2.42) 
for f3 7^ a and 

»K(*a,r a ,o;) = ^Pi, 2 (r a ,a a ,a;)(n a (a;) x <J, (2.43) 

»2^(ila,r a ,a;) = £ T^(r a ,a a ,^) • / a ,, 2m », (2.44) 

PSU^'a,") = £ DS(r a ,a Q ) • / Q ,, 2m » (2.45) 

for (5 — a. 
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Here, R a p = R a — Rp = (R a p, ®ap, ®ap) is the vector between the centers of spheres a and 
(5 pointing from sphere j3 to sphere a in the spherical coordinate system O. The quantities 
f/3,im( u ) are the harmonics of the induced surface force densities fp{ap,uj) expanded in the 
spherical harmonics 

1 



(2.46) 



fp{a p ,u) = — =J2f pM (ap : u;)Y lm (9p : (pp) : 

V lm 

where, for short, the argument ap in fp i m (ap, oS) is omitted. The quantities involved in 
Eqs. (2.40)-(2.45) are written as 



lm 



l 2 m 2 



2 f k 2 

F h i 2 ,i{r a ,a l3 ,R a p,u) = — Jdk 



k 2 + k- 



■ jh( kr a)ji 2 (kap)ji(kR a p), 



limi,lm 



h-h+l 



JdVL k (I - n k n k ) Y* imi (0 fc , (p k )Y hm2 (9 k: (p k )Y* m (9 k , tp k ), 
D\lZ(ra, ap, R a p) = £ *>fe m (r«, ap, R a p)Y lm (&«p, 

lm 

oo 

2 r 

C h h,i( r c*, a p, R ap) = - dkk j h (kr a )ji 2 (ka )ji(kR a(3 ), 

7T J 

1 yiifi fc n fc Y* imi (9 k , <Pk)Y hma (9 k , <Pk)Y^(9 k , tp k ), 

oo 



W l ? m2 , - 1 

l\m\,lm 





h-h+l- 



2 r k 

Phh,i( r <*i a pi R ap,v) = — jdk 2 j h (kr a )j h (kap)ji (kR a p) . 



T l2 Zl(r a ,a a ,u; 
F h i 2 (r a , a a ,u 



^hml( r a-> a a 



Phh(r a ,a a ,u 



1 



-^j^ F h i 2 (r a , a a , v)K\l™ 2 uQO i 



2 f k 2 

ITT) J k z + K 
' 



1 



oo 

2 r 

= - dkkj h (kr a )ji 2 (ka a ), 

7T J 



2 f k 

= dk 71~, ojh( kr a)jl 2 ( ka a)- 

TYTj J k 2 + K 



2.47) 
2.48) 
2.49) 

2.50) 

2.51) 
2.52) 
2.53) 

2.54) 

2.55) 

2.56) 
2.57) 
2.58) 
2.59) 

2.60) 
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The explicit form for the quantities defined by (2.49), (2.50), (2.52)-(2.55), and (2.57)- 
(2.60) has been found early in [39] (see Appendix A). 

The quantities e m , where m — 0, ±1, are defined as follows: 

e = e z , e±i = (ie y ± e x ) , (2.61) 

i.e., up to the factor (— l) m , are the cyclic covariant unit vectors [41], and e x , e y , and e z are 
the Cartesian unit vectors. 

Since the quantities K\ 2 ^ i lm and lm are not equal to zero, respectively, only for 

I = h + l 2 -2p, where p = -1,0, 1, . . .,p max , Pmax = min([(/i + l 2 )/2] , 1 + min(/i, / 2 )) , where 
[a] is the integer part of a and min(a, b) means the smallest quantity of a and b, and for 
l = h + l 2 -2p+l, where p = 0, 1, . . . ,p ma x, Pmax = min([(Zi + l 2 + l)/2] , 1 + min(Zi, Z 2 )), 
the infinite sums over / in relations (2.47) and (2.51) are replaced by sums containing a finite 
number of terms corresponding to these values of I. The quantity 00 ^ only for 

l 2 = h ± 1 > 0. Therefore, the infinite sums over l 2 in relation (2.40) and (2.45) are replaced 
by sums containing only terms with l 2 — l\ ± 1 > 0. Taking into account that -K^^ oo 7^ 
only for l 2 = h + 2p > 0, where p = 0, ±1, the infinite sum over Z 2 in relation (2.44) is 
replaced by the sum with l 2 = li + 2p > 0. 

According to (2.43), only one harmonic of v^i^(Ra,r a ,u) is not equal to zero. Using 
this observation, we can represent this component of the fluid velocity as follows: 

v^ md (R a + r a ,uj) = ^P li2 (r a ,a a ,u)(p a (u) x n a ), (2.62) 

where n a = r a /r a . 

Relations (2.29), (2.31)-(2.33), (2.36), (2.37), (2.40), (2.41), (2.43), and (2.44) completely 
determine the function v(r,u) in the entire space provided that the induced surface force 
densities are known. Analogously, relations (2.30), (2.34), (2.38), (2.42), and (2.45) uniquely 
reproduce the original function p(r, u) at any point of the space with the exception of 
the points of the surfaces of the spheres (r Q = a a , a = 1,2, ...,N) where the obtained 
quantity is equal to the half-sum of the original function p(r,cu) given at r — R a + a a + 
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and r = R a + a a — because the original function p(r,co) defined by relation (2.13) is 
discontinuous at the surfaces r a = a a . The fluid pressure at the surfaces of the spheres 
denoted as p(R a + a a + 0, a;) can be represented as follows: 

p{R a + a a + 0,u)= P {0) {R a + a a , u) + p ind (R a + a a + 0,u), (2.63) 

where 

p ind (R a + a a + 0,cu) = 2p ind (R a + a a , u) - icop (R a + a a ) • U a (u) (2.64) 

and the quantities p^(R a + a a ,u) and p md (R a + a a ,u) are defined by the corresponding 
expressions for p^(r,u) and p md (r,u) given at r = R a + a a derived above. 

Using the expansion (2.46) together with Eqs. (2.23) and (2.24), one can represent the 
Fourier transforms of the force F a (u) and the torque T a (uj) exerted by the fluid on sphere 
a in terms of Fourier harmonics of the induced surface force density f a i m (w) as follows 
[34,35,39,42]: 

f-M = -/«,ooH, (2-65) 
T» = E ( e ™ x /«,imM) • (2.66) 

V<J m =-l 

Thus, both the fluid velocity and pressure induced by the spheres as well as the forces 
and torques exerted by the fluid on the spheres are expressed in terms of the harmonics of 
the induced surface force densities. Note that these expressions have been obtained without 
imposing any additional restrictions on the size of spheres, distances between them, and the 
frequency range. To determine these harmonics, using the stick boundary conditions for the 
fluid velocity at the surfaces of the spheres [8,13,14], the original problem is reduced to the 
determination of the unknown quantities fp i m {oS) from the following infinite system of the 
linear algebraic equations [39]: 

E T a a tZ» ■ fa,i 2 mM + <S2» - v a , hm M = - E { E t^z\{u) ■ f PMm M 

+ <&>)}, ( 2 - 67 ) 
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where 

V aMmi (u) = 5i lfi 5 mifi U a (uj) + 5 hA (n a (w) x e^J , (2.68) 
and, to simplify the representation, we introduce the following notation: 

where the quantities T^(a Q , a^, i? a /?, w), T^(a a , a a , w), u^j^(i2 a , a Q , u;), and 
v c*2.mi (-^aj a «, w), are defined, respectively, by relations (2.47), (2.56), (2.40), and (2.43) 
at r a = a a . Up to the terms v^i^(R a ,a a ,u), where (3 = 1,2,...,N, Eqs. (2.67) agree 
with the corresponding equations given in [34,35] (for the discussion, see [39]). Follow- 
ing [34,35,42], the quantities T^^aa^p, R a p,u), where f3 ^ a, and T l fi™l(a a ,a a ,iv), are 
called hydrodynamic interaction tensors (mutual interaction tensors for (3 ^ a and self- 
interaction tensors for (3 — a). The quantities Fi^ifaa, ap, R a p, u) and Fj^^aa, a a , u>) defin- 
ing the hydrodynamic interaction tensors, the quantities Cy 2 ,/( a co ap, R a p) and C^j^aa, a a ) 
defining, respectively, the tensors D l l ^^(a a , ap, R a p) an d £)^™J(a a , a a ), and the quantities 
Pi l2 ,i(a a ,ap,R a p,u) and Pi, 2 (a a , a a , a;) defining, respectively, v^ mi (u) and ^^(w) are 
determined by relations (A6), (A16), (A7), (A10), (A14), and (A18) at r a = a a . For (3 ^ a, 
these quantities have the form 



Fi^iia^a^Ra^u) = — { jh(K)ji 2 (bp)hi(y a p) - 5 lM+la+2 



Clih,l\ a <xi a /3i Rap) — $l,h+li+l 



^ I 2 r (h + 1) r (i 2 

< i 

r(z 1 + z 2 + 1) 



x _^Ja^ / = / 1 + / 2 _ 2 p>0, p = -l,0,l,...,?w, (2.69) 



2 r(i 1 + f)r(z 2 + §) Kp ' 

Z = Zi + l 2 -2p+ 1 > 0, p = 0,l,...,p 



maxi 



P h2 Ma,ap,R a p,Lu) = (-l) p —jiMUh)Wy a t)), l = h±l>0. (2.70) 

711] 



Here, j;(x) = yir/(2x)I l+ i(x) and = y n/(2x)K l+ i (x) are the modified spherical 

Bessel functions of the first and third kind, respectively, [43], T(z) is the gamma function, 
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the dimensionless parameters a a p = a a j R a p and ap a = ap/R a p, which are the ratios of the 
radius of a sphere to the distance between the centers of two spheres, is always smaller than 
one (for spheres of equal radii, they are equal and cannot be greater than 1/2; in dilute 
suspensions, cr a p,ap a <C 1), the dimensionless parameter y a p = nR a p (\y a p\ = V2R a p/5) 
characterizes the ratio of the distance between the centers of spheres a and (3 to the depth of 
penetration 5 of a plane transverse wave of frequency cu into the fluid. In certain frequency 
ranges, the quantity \y a p\ can be both smaller (for short distances between spheres) and 
greater (for space-apart spheres) than one. 
For p ^ a, we have 

F hh (a a ,a a ,u) = (-iy^j lmax) (b a )h lmin (b a ), l 2 = h + 2p>0, p = 0, ±1, (2.71) 
C h i 2 (a a ,a a ) = — (Si 2ih+1 + 5j 2)Jl _i) , k = h ± 1 > 0, (2.72) 

2 ~ 

Pi,2(a a , a a , u) = — h 1 (b a )j 2 (b a ), (2.73) 

717] 

where l max = max(/i,/ 2 ) and l rnin = min(/i,/ 2 ). 

The agreement between the quantities F^^aa, ap, R a p,uS) and Fi li2 (a a , a a , uS) defined 
by relations (2.69) and (2.71) and the corresponding results given in [35] is discussed in 
[39]. Here, we only note that, in the general case, according to (2.71), the self-interaction 
hydrodynamic tensors are not equal to zero for three values of l 2 [l 2 — h, Zi + 2, 

and Zi — 2 (for l x > 2)] instead of two values of l 2 [l 2 = h, h + 2] mentioned in [35]. 

In [34,35], the similar system of equations in the unknown harmonics of the induced 
surface force densities is considered for equal spheres (a a = a, a = 1, 2, . . . , N) for 
b a = b, \b\ 1 and the corresponding solution of the system is given as a series in two di- 
mensionless parameters (a and b). However, the expressions for the harmonics of the induced 
surface force densities as well as for the forces and torques exerted by the fluid on the spheres 
determined up to the third order of these parameters are expressed in [34,35] in terms of 
the tensors T„'[^(w) inverse to the self-interaction hydrodynamic tensors T„'{^(w) and 
defined by the condition 
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£ CM'CM^'.*^,^' (2-74) 



m 3 =-« : 



In [35], several inverse tensors [e.g., T^^iu) and T„^j(w)] are expressed in terms of 
the inverse tensor K lm 2 i 00 , where the tensors X;| m ^ 00 inverse to the tensors ^^,00 are 
defined by the condition 

£ ^imi,00 ' K-llml, 00 = ^i.mj^i /i = 0, 1, 2, . . . . (2-75) 

m3=— h 

However, the explicit form for the tensors -K^^oo ^ s 110 1 gi ven except for li = 0. At the 
same time, as is shown in [39], 

detK^ j00 = 0, (2.76) 

which means that the tensor 00 is singular, and hence the inverse tensor -K~j™^ 00 does 

not exist (at least, in the conventional sense). As is shown below, this difficulty remains 
valid even in the case of one sphere. 

In the present paper, we use the procedure proposed in [39] which is free of the above- 
mentioned difficulties. According to this procedure, in addition to the original system of 
equations (2.67), the following equations may be used: 

Y a (Lu) + iup4nV3al(R a -U a (iu)) =4irV3alY,P l p d (R*,u), a = 1, 2, . . . , TV, (2.77) 
where 

Y a {w)= £ e m -/ Q)lm H, a = 1,2,..., N, (2.78) 

m=— 1 

i /• 00/2 ?2+l 

pf\R a ,u) = — 5- { (n^ • /aooM) + 4tt £ £ £ w*S£ +1>m . / AIama M 

47r/t a/3 k / 2 =lm 2 =-/2 m=-(Z 2 +l) 

x^ 2+ i iro (e Q/ 3,$ a/3 )}, (2.79) 
and n a/ 3 = R Q p/R a p. 

For the derivation of Eqs. (2.77), (2.79), it is necessary to set r a = e, where e — > +0, in 
relations (2.30), (2.34), (2.38), (2.42), and (2.45), take into account the explicit form for the 
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function p(r, u;) denned by relation (2.13) and relations (A19) and (A20) for the quantities 
Qi« 2 a 0-> Rad) an d lim C; 1 ; 2 (r Q ,, a a ) given in the Appendix A. 

We investigate the system under consideration assuming that \b a \ <C 1 (a a <C S/y/2 , i.e., 
\uj\ <C vjdfy and represent all required quantities [the harmonics f p t i 2Tn2 {oj) of the induced 
surface densities, the fluid velocity and pressure, the forces and torques exerted by the fluid 
on the spheres, and the friction and mobility tensors] as power series in two dimensional 
parameters (a and b) retaining terms up to the order of b 3 . 

In this approximation, the quantities F^^Oo,, a^, i? aj g, u;) and Fj 1 j 2 (a a , a a , uj) defining 
the hydrodynamic interaction tensors are simplified to the form 



rjR a p r (h + 1) r (l 2 + 1) 1^2^^ 

i ( K b% 

32 \{h + §)(«! + 1) ('2 + |)(i2 + §) 



2 



1 , U bl , b 



4V/i + | / 2 + 3 
- ^t! /. ,w. ^ ,w. -,, ) • Oil!', 



hi{y a p) 



r (h + h + z . , 

- 2 - L \,l = li + k-2p,p= -1,0,1,..., Pmax , (2.80) 

Va/3 ) 

F hl2 (a a ,a a ,uj) = F h (a a )P h>n (b a ), l 2 = l 1 + 2n>0, n = 0,±l, (2.81) 
where the quantity F^(a a ) corresponding to the stationary case is defined as follows: 

F '-(-»-(^TT)^ (2 ' 82) 



and 



Pi,M = l-b a d l (b a ), (2.83) 



, 1 — 5;,o b a 



p '-^ - -"'w^m^-y (2 - 85) 

d,(b a ) = 6,,„ + b a _ 2 + + 1 (i,, - S hl ), (>0. (2.86) 

We expand the quantities Pi,2(a a , a a , uj) and -Pz^z^a, R a /3,^) in power series in the 
small parameter b retaining terms up to b 3 and substitute these expansions into relations 
(2.43) and (2.40) at r a = a a . As a results, we get 
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a n b 



<K» = *l„i ^ x e;„,) + (2.87) 



1 ^i + ^J l 2 =d±l>l m 2 =-l2 



x Mva^i&^e^ , + o{b A ). (2.88) 

For the solution of the system of equations (2.67), (2.77), we take into account relations 
(2.79), (2.87), and (2.88) and the following estimates for the hydrodynamic interaction 
tensors following from (2.80) and (2.81) valid for any l±,l 2 > 0: 

a a T%lZl M ~ <T° aP i>a, h h (ba), h = h + 2n > 0, 71 = 0, ±1, 

aaT^Zl M ~ rtp^X^Mh (&«, bp, y aP ) , (2.89) 

where ip a ,hh(ba) an d ipapjihibonbpjyap) are certain finite functions of the small parameters 
b a and bp and the parameter y a p free of any restrictions imposing on it. The function 
^a,hi 2 {ba) satisfies the relation 

\imJ ahl2 {b a ) = 5 l2M a a T%\Z\, (2-90) 

where 

rpa,hm,2 V^" jxlim2 /o qi \ 



is the static hydrodynamic self-interaction tensor [35,39,42]. 

These estimates enable us to seek a solution of system (2.67), (2.77) by the method of 
successive approximations using Eqs. (2.67), (2.77) with the zero right-hand sides as the 
zero iteration, which corresponds to the absence of interactions between the spheres. This 
is quite natural because the sum and the second term on the right-hand sides of Eqs. (2.67) 
correspond to the harmonics of the fluid velocities at the surface of sphere a induced by 
the surface and volume force densities distributed, respectively, over the surface of sphere (3 
and inside its volume, whereas analogous quantities on the left-hand sides of these equations 
correspond to the harmonics of the fluid velocities at the surface of sphere a induced by the 
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surface and volume force densities distributed, respectively, over the surface of this sphere 
and inside its volume. Analogous reasonings may be used for Eqs. (2.77). 
Next, we represent the unknown harmonics f a j m (uj) in the form 

oo 

/ q mH = £/S>), (2-92) 

n=0 

where f^a) imi (ou) ls a solution of system (2.67), (2.77) corresponding the nth iteration. As 
is shown in what follows, these quantities can be represented in the form 

/S» = + /SiM + hi filial n = 0, 1, 2, . . . , (2.93) 

where fafmi 1 ^) anc ^ fa'^li 00 ) are the components of the harmonic f^Jmi 10 ) associated, re- 
spectively, with the translational motion of sphere a and its rotation and f^hli 00 ) i s the 
harmonic of the potential component of the induced surface force density. In what follows, 
we show that the last quantity has no effect on the fluid velocity and the forces and torques 
exerted by the fluid on the spheres. 

In the theory of hydrodynamic interactions between spheres, the main problem is to 
derive the required quantities up to the terms of a given order in certain small parameters 
(in the general nonstationary case, these are the sets of parameters a and b). Using the 
estimates (2.89), it can be shown that, for the nth iteration, the main terms of the power 
expansions of the quantities fa'jmi 00 ) an d f^afm( u ) * n a are proportional to a l+n and a l+n+1 , 
respectively. This enables us to formulate the following procedure for determination of the 
induced velocity and pressure of the fluid (as well well the forces and torques exerted by the 
fluid on the spheres) up to terms proportional to a p , where p is a certain positive integer: 

(i) to carry out p iterations; 

(ii) for each sth iteration, s < p, to keep only harmonics with I — 0,1, ... ,p — s; 

(iii) for each keeping harmonic, to retain all terms of power series up to terms proportional 
to a p inclusive. 

For the fluid velocity and pressure, further simplification of the approximate results is pos- 
sible depending on the choice of a point of observation. 
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In the present paper, we restrict ourselves to the consideration of only the first two 
iterations. For this, it is sufficient to retain the terms of at most the second [for fj£\w)] and 
the third [for f^\u)} orders in a, respectively. 

Using Eqs. (2.8), (2.9), (2.21), (2.22), (2.27), (2.65), (2.66), (2.92), and (2.93), we repre- 
sent the total force F^fa) and torque T^°*(c<j) acting on sphere a as follows: 

F tot {u) = _ F* {u) + {2M) 

T tot {u) = T e : t {u) _ r ,t {u) + To(w)j (2 ^ Q5) 

where 

^ a (u) = Ff(u) + Fi\u) 1 (2.96) 

T a (u) = T®{w) + T^\u) (2.97) 

are the force and torque exerted by the fluid on sphere a due to the translational motion 
and rotation of all spheres, 



^ ) M = n n M + E^'" ) M, (2.98) 

F™{uj) = -ium a U a {uj), (2.99) 

oo 

*l r) M = E*t ,n) M, (2.100) 

n=0 

oo 

^ ) M = E^' n) M, (2.101) 

n=0 

oo 

T a\u) = 5Z T i r ' n) H> (2-102) 

n=0 

and F^ n) (uj), F^ n \oj), T^ n) (uj), and T^ n \uj) are the forces and torques caused by the 
translational motion [superscripts (t, n)] and rotation [superscripts (r, n)] of the spheres 
corresponding to the nth iteration defined as follows: 

*t C ' n) M = -/&o } H, C = t,r, (2.103) 

Tp } H = ~ E (c™ x , C = t,r. (2.104) 

V<J m=-l 

The similar representations take place also for the induced velocity and pressure of the 
fluid. 
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III. NONINTERACTING SPHERES 



In the approximation of noninteracting spheres (zero iteration), the infinite system of 
equations (2.67) is splitted into the collection of independent systems of equations for each 
sphere a, 




a = 1,2,..., TV. (3.1) 

Note that, for each a, we still have the infinite system of linear algebraic equations in the 
unknown quantities faj m (u)) (instead of a collection of systems for each value of I in the 
stationary case [39]). 

The system (3.1) may be solved together with Eqs. (2.77), which are rewritten in the 
approximation of noninteracting spheres as 

Y^\uo) = -iup4nV3al(R a -U a (u)), a = 1, 2, . . . , N, (3.2) 

where 

Yi°\u>)= £ e m -/i 0) lr >), a = 1,2,..., TV. (3.3) 

ra=— 1 

Physically, the presence of the non-zero term on the right-hand side of this equation cor- 
responds to the appearance of the potential component of the induced surface forces even 
in the approximation of noninteracting spheres (at the expense of an account for the time 
dependence of the problem). 

In the general case of interacting spheres, the conditions of the absence of the potential 
components of the induced surface forces have the form 

Y a (u) = 0, a = 1,2,..., TV. (3.4) 

Solving the systems of equations (3.1) and (3.2) up to terms of b^, we get 



23 



= ft,o /£S!m + /Sm, (3.5) 

2a r 



= <V /SH = A>)(n a M x <), (3.6) 



= <4^A m = 0, ±1, (3.7) 



where 



/JolM^^fM^M, (3.8) 
= ^K m .[/ a ( W ), m = 0, ±1, ±2, (3.9) 



p(b) = l + b + j, (3.10) 

cM = o r (y, (3.H) 

being the Stokes stationary friction coefficient for a sphere rotating in a viscous fluid 
(C a = S^al), and 

s r (6 a ) = l + ^(l-6 Q ) (3.12) 

is the corresponding frequency factor. In deriving (3.9), the use of the tensor K m = 

2V307T i^2m,oo> where m = 0, ±1, ±2 (see Appendix A), and the relations 

2 _ 2m 

X! ^2mi ,00 ' ^2m 2 , 00 = -^2mi,00 (3.13) 

m2=— 2 

have been made. 

Note that in the approximation of noninteracting spheres, the expression for the harmonic 
fa'Jmi ^) gi ven i n [35] is proportional to the tensor ^T lm ^ 00 that is inverse to the singular 
tensor K\™* 00 , whereas according to (3.6), the harmonic fa'iL( u ) contains no singularities. 

Substituting the time-dependent solution (3.5)-(3.7) into Eqs. (2.41), (2.42), (2.44), and 
(2.45) and using (3.13), we obtain the following relations for the harmonics of the fluid 
velocity and pressure induced by the surface forces in the approximation of noninteracting 
spheres: 

V aMm7 d ( R <*i r <*i UJ ) = —P( b <*) I <^i,0 <W,0 3o(b a ) exp(-X a )U a (w) 



J a,hm 1 

1 a 
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+ S h , 2 - - 3^ h 2 {x a )Ub*))K mi ■ U a {u)y (3.14) 

v *hu}r d ( R ^ r ^ UJ ) = Ki^^C a (u)Fi,i{r a ,a a ,uj) (p a (w) x e^J , (3.15) 
vff^; nd ( R a,r a ,u;) = ^p(b p )T™ mi (r a ,ap,R aP ,u;) ■ Up(u), (3 ± a, (3.16) 
vflT^ 11 *, r a , *) = -?- CM E E ^i,/(r Q , oa, i^, a;) (O^a;) x Wg^) 

ZO /3 J=h±l>l m=-l 

x Yi m (e a/J ,$ a/J ), /3^a, (3.17) 

«Sfm! ind («a> r a , w) = 0, P = 1,2,...,N, (3. 18) 

Pi^l^^a^) = " £ (l + &« + £) (e^ 1 ■ £/»), (3.19) 

PpSZ d ( R <*> r «> w ) = &P(h) D h mi ( r «> ^) • Up{w), (3^ a, (3.20) 
P^A^- r a ,u>) = 0, (3 = 1,2,-.., N, (3.21) 
r a , u) = 5 p , a S lu0 5 raAa % - up (R a ■ U a {u)) , f3 = 1, 2, . . . , N. (3.22) 

Note that, by virtue of (3.21), the rotation of the spheres has no effect on the fluid 
pressure, 

P $ fi * nd (r, U ) = 0, (3 = 1,2,...,N. (3.23) 

In view of (3.7), the potential component f^'°\a a , u) of the induced surface force density 
f£\a a ,u) in the approximation of noninteracting spheres has the form 

f r\a a ,u) = ^^M, (3.24) 
a a 47tv3 

and, hence, has only the radial component. As a result, the potential component of the 
induced surface force has no effect on the forces and torques exerted by the fluid on the 
spheres immersed in it. According to (3.18), the potential component f^'°\ap,u), where 
(3 = 1,2, ...N, of the induced surface force density makes no contribution to the fluid 
velocity. 

These conclusions concerning the influence of the potential component of the induced 
surface force densities on the fluid velocity and on the forces and torques exerted by the fluid 
on spheres also remain valid for any nth iteration if the potential component f^' n \a a ,u) 
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of the induced surface force density f^'(a a ,uj) corresponding to this iteration has only 
harmonics with I = 1 defined by a relation similar to (3.7). 

According to relations (3.19) and (3.20), for the fluid pressure for r a > a a , the account of 
the time dependence leads only to the appearance of the frequency-dependent factors p(bp) 
equal to 1 in the stationary case and the substitution of U p(uo) for Up. 

Taking relations (2.87) and (3.15) into account, we obtain the following expression for 
the fluid velocity: 

v^ nd (R a + r a , u) = vf^ m \R a + r«, u) + v^ ind (R a + r a , u) 

= (—) q(K)ip(xa) (n a (w) x r a ), (3.25) 

where 

b 2 b 3 

9(*) = l + y-J> ( 3 - 26 ) 

^(x) = -x 2 h 1 (x), (3.27) 

7T 

at an arbitrary point r = R a + r a induced due to the rotation of sphere a, which agrees 
with well-known result [14] for the fluid velocity induced by a rotating sphere for \b a \ <C I. 

Using relation (3.19), we obtain the following expression for the fluid pressure for r a > a a 
caused by the translational motion of sphere a: 

p (t,o)ind {Ra + rajUj) = Z JL^ p (p a ) U . U a (u)\ (3.28) 

which differs from the corresponding relation for the stationary case [14] only by the 
frequency-dependent factor p(b a ). 

If the point of observation is far from sphere a (Rex a ^> 1), then using relations (3.14) 
and (3.25) and retaining the leading terms, we obtain the following expressions for the fluid 
velocity induced by the motion and rotation of sphere a: 

v {Smind {Ra + ) w 3o«p(M ^ _ 

2r a x 2 a 

\T a J 
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Thus, the account of the time dependence in the Navier-Stokes equation leads to the essential 
change in the asymptotics of the fluid velocity in the far zone [14,16]. [Of course, here and 
below, it is assumed that r a is such that the linearized Navier-Stokes equation (2.1) can be 
used. For details, see Appendix B.] 

Next, using relations (3.16) and (3.20), we obtain the following expressions for the fluid 
velocity and pressure induced in the far zone (r Q ^> R a p) due to the motion of sphere (3 ^ a 
with the velocity Up(oj) in the zero approximation with respect to the small ratio R a /3/r a - 
v (s,t,o)in d{Ra + ^ w) w Z^P(M_ {3nana _ J} . (3 31) 

pf 0)ind (R a + r a ,u>) « ^ %p{b p ) (n a ■ U p {u)) + 0(0%). (3.32) 

The use of Eqs. (3.28), (3.29), (3.31), and (3.32) results in the following space distribu- 
tions of the fluid velocity and pressure induced by the motion of iV spheres far from them 
in the zero approximation with respect to the ratio R a p/r a : 

o N 

v (s,t,o)ind {Ra + r ^ u) w _ ( 3riana _ j) . £ a p p(bp)Up(u), (3.33) 

o N 

p {t,o)ind {Ra + rajUj) ^V^ apv{bp) (n a . U p (u)). (3.34) 

Zr a 13=1 

In the case of equal spheres ap = a (bp = b), these relations are simplified to the form 

v (S,t,0)ind( Ra + u) „ 3_ Qp(6) (g^^ _ j) . (3>35) 

p (t,o)ind( Ra + ^ w | ^ ap ( 6 ) ( nQ . E/ tot (u;)), (3.36) 

where 

U to \uj) = Y J U^). (3.37) 

0=1 

Hence, the velocity and pressure fields of the fluid caused by the translational motion 
of equal spheres in the far zone, in the zero approximation with respect to the ratio of the 
typical distance between two spheres to the distance to the point of observation, can be 
represented as the velocity and pressure fields of the fluid induced by a single sphere of 
radius a moving in the fluid with the translational velocity U tot {t). 
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If all spheres move with the same translational velocity Up(u) = Uo(u), we get 



a 



N 

where a e tt = J2 appibp). Thus, in the far zone, the action of the system of spheres moving 
with the equal velocities is equivalent to the action of a single sphere of radius a e ff. 

Using (2.103), (3.5), (3.6), (3.8), and (3.9), in the approximation of noninteracting 
spheres, we obtain the following relations for the Fourier components of the force exerted 
by the fluid on sphere a caused by translational and rotational motions of the spheres: 

F^\u) = -Zap(b a )U a (u>), (3.40) 
*t 0) M = 0. (3.41) 

If we restrict our consideration to the approximation of noninteracting spheres, then, 
using relations (2.98), (2.99), and (3.40), we obtain the following expression for the force 
acting by the fluid on sphere a moving in it with the velocity U a {uj) 

?M = -U")U a (w), (3.42) 

where 

U") = ^s\b a ) (3.43) 

is the Stokes frequency-dependent friction coefficient for a sphere of radius a a moving with 
the velocity U a (u) and 

s\b a ) = l + b a + h ^ (3.44) 

is the frequency factor corresponding to the translational motion. 

Equation (3.42) defines nothing but the celebrated Boussinesq formula [14,16] for the 
force exerted by the fluid being at rest at infinity on a single sphere moving with time- 
dependent velocity. Indeed, returning to the time variable, we immediately obtain 

28 



<zr u\ c tt u\ ™*dU a (t) 2 f dr dU a (r) 

— oo 

To determine the spectral component of the torque exerted by the fluid on sphere a, we 
use relations (2.104), (3.5), and (3.6) and obtain the well-known results [31] 

T%°\u) = 0, (3.46) 

rf' 0) H = (3.47) 

In the time representation, this reads 

T (t) - -m (t) - 5/ dQa{t) - / dT d " a(r) (3 48) 

— oo 

where J Q = (2/5)m a a^ is the moment of inertia of a fluid sphere of radius a a . Analogously 
to the Boussinesq force (3.45), torque (3.48) exerted by the fluid on a sphere depends on 
the angular velocity of the sphere, its acceleration, and the history of rotation of the sphere. 
Furthermore, the second term in (3.48) is independent of the fluid viscosity. 

Finally, using the relations (2.5), (2.6), (2.94), and (2.95), we can represent the equations 
of motion of sphere a as follows: 

[-iurrgf + £ a (l + b a )}U a {u) = Ff( W ) - F e a x \u), (3.49) 
[-iu{l a + 5(1 - b a )I a } + £}n«(") = TfH - t;f H, (3.50) 

where rrfj* = m a + m a /2 is the effective mass of sphere a [14,16]. 

IV. THE FIRST ITERATION 

For the first iteration, the system of equations (2.67), (2.77) is reduced to the form 

E ■ /£L» = - E { E rJSM • /gL» + ( 4 - x ) 

YW(u>)= j: e m .f^ m (uj) = AnV3alY,pf md (R a ^), " = 1,2,..., A, (4.2) 

m=—l p^a 
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where (R a ,ui) is defined by relation (2.79) with fpj 2m2 

■^ftl^i^)' which gives 

y«M = VS E • E/^H), a = 1, 2, . . . , N. (4.3) 

As before, the necessity in Eq. (4.3) is caused by (2.76). Solving Eqs. (4.1) and (4.2) 
up to the terms of the order b 3 and a 2 [for /^r^O^)] an d c 3 [for /^ m2 ( a; )] anc ^ usm § 
Eq. (3.13), we obtain 



/££M = *i,o /So'H + <V + *i, 2 /g&M + *, >3 /SglH, C = f , r, (4.4) 

= < m = 0, ±1. (4.5) 



Here, 



/SoM = -£«p(&«) E tpp{bp)T M {R aP , u) ■ Up(cu) + 0(<x 3 ), (4.6) 

/¥<* 

/£$(<") = E apofcqibpMyafdfnafl x fyH) + 0(<r 4 ), (4.7) 



/£& = -f & + (6^H) T } , C = t,r, m = 0, ±1, (4.8) 

= • E(! +M^m(^,o;) • t^w), m = 0,±1,±2, (4.9) 

V3U ^ a 

= m E 3 ^ • {4c<S» + (c£S>)) T }, 

m = 0,±1,±2,±3, (4.11) 

m = 0,±1,±2,±3, (4.12) 

*£2M = *££ 0) M + + *2#V), ™ = 0, ±1, (4.13) 

&£#M = &S 0) H + ^(u;), m = 0, ±1, (4.14) 

&g£ 0) M = E &S 1,0) M, C = t,r, m = 0,±1, (4.15) 

^ 1,0) H = ( W ) • U p (u), m = 0, ±1, (4.16) 

&&!' 0) M = 2^ E E ^K,a^,it:^,a;)(^H x ^J^fe^^^), 

1=0,2 m=-l V 7 
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mi = 0, ±1, 



(4.17) 



+ ( b <t,i,o), 



CO 



e ) M = E^ 1 ' 0) M, m = 0,±l, 



1) 



771 1 

m = 0,±1, 

T 



(4.18) 

m = 0,±l, (4.19) 
(4.20) 



+ 2& E E Pi, 2> i(aa,a/9,i^,a;)(n j9 (a;) x < 1|00 1^(6^, $ Q/3 ) 



+ 6 



Z=0,2 m=-Z 

2 4v^F 



a 15 



(i-jk) £ k 

77l2 = -l 



17T12 
17771,00 



46S 2 °H-) + (6S 2 )H) T }, 
mi = 0, ±1, 



c£SM = &S 0) M + I<g(a;) = E( X + M^H' = 0, ±1, 



b 2 b 3 

r(b) = l + b + - + -, 



(4.21) 
(4.22) 

(4.23) 



+ 



2 ^o,o,2(a a , a/3, -Rq/3, (3n a/ 3n a/ g — l) | 



(4.24) 



is the modified dynamic Oseen tensor [34,35,39]. This tensor is symmetric both in the 
coordinate indices = x,y,z) and the indices of spheres (a,/3). Note that in the case 
under consideration, only terms corresponding to the approximation of the problem with 
respect to the small parameters b and a should be retained in Tm(Ro>/3, w), namely, 

1 



-^0,0,0 (d a , Up, Ra/3, 
^0,0,2 (Go:, dp, Ra/3, w) 



r]R af 3 

1 



1 + a ex P(-2M), 



(4.25) 



a/3 



1 + 



bl + b} , b% + b% 



6 



120 



2 q r 



-irMvafi) ■ (4-26) 



7T 



In deriving relations (4.11) and (4.12), after the determination of the inverse tensor 
lC,oo, the tensor 



31 



N mim2 =2J 7 - ]T k 3 Zl 00 ■ KlZloo, mi = 0,±l,±2,±3, m 2 = 0,±l, (4.27) 

V <J m3 =-3 

has been introduced. Its components are expressed in terms of the tensor K m with m = 
0, ±1, ±2 given by the relations (A2) in Appendix A as follows: 



N mm — |5 m ,o + (l — S mfi ~^K , 



m = 0,±1, 



2 1 1 

i>m,0=^Am, "l = ±l, JV , m = -j= K _ m , 171 = ±1, iV m _ m = -j= K ± 2 , 771 = ±1, 

N ±2fi = ^K ±2 , N 2il = ^K u N^^sj^K^ N 3A = ] j^K 2 , 

JV-3,-1 = ]J^K-2, JV 2 _i = AT_ 2)1 = iV ±3 , = iV 3 ,_i = JV_ 3) i = 0. (4.28) 

Here and below, the three superscripts in quantities b^'l^\uo), where ( = t,r and n,k = 
0, 1, 2, ... , mean the following: £ shows that the induced surface force densities are caused 
by the translational ( = t) or rotational (£ = r) motion of the spheres, and the subscripts 
n and fc stand for the corresponding iterations with respect to the parameters a and b, 
respectively. 

All quantities of the type (b a ^ m ) T are defined as follows: Consider three vectors b m = 
(b mtX ,b mj y,b mtZ ), where m = 0, ±1, and decompose them into the independent unit vectors 
e and e±i defined by relations (2.61) taking into account that any vector a = (a x ,a y ,a z ) 
may be represented in terms of these unit vectors in the form a = (a +1 , a_i, a ), i.e., 

a = — a-\e\ — a + ie_i + a e , (4.29) 

where 

a ± i = ^|^. (4.30) 

Accordingly, we introduce the vectors 6 m = (6 m +i, 6 m -i, 6 m o) and represent them as 

i 

& mi = S mi m 2 -e m2 , TOi = 0,±l, (4.31) 

m2=— 1 

where 

32 



B 



boo 


-b -i 


-b +i 


b+w 


-b+i-i 


-6+1+1 


b-w 


-6-1-1 


-6-1+1 



V 

Then is the vector defined as follows: 

i 

mi — - D mirrt2 ' e ™2 J 

mi=— 1 



mi, m 2 = 0, ±1. 



(4.32) 



mi = 0, ±1, 



(4.33) 



where Bj nim2 is the matrix transposed to the matrix B mim2 . 

Substituting (4.6) and (4.7) into (2.103) and (4.8) into (2.104), we obtain the following 
relations for the forces and torques for the first iteration: 









(4.34) 








(4.35) 








(4.36) 


r^(a;) = 


-Etf'M 




(4.37) 



where £^3^ (a;) and (a;) are, respectively, the translational and rotational friction tensors 
corresponding to the first iteration and an d ^2^\ w ) are the friction tensors that 

couple translational motion of spheres and their rotation. These quantities have the form 



= -U0(b«)p(bp)T M (R a p,uj) + 0(a% 
t^iPiv) = -^pa j3 a a i3ap a T{b Q )q{bp)il}{y ol p) (e • n af3 ) + 0(a A ), 
i r S\ w ) = -^.a^ a i3(ypaT{b j 3)q{b a )'4){y al3 ) (e • n a/3 ) + C(cr 4 ), 
Cf\co) = ^a 3 a q(b a )q(b p )y 3 ap {h 2 (y^)(l - 3n a/3 n Q/3 ) + 2h (y a p)l} + G(a 5 



(4.38) 
(4.39) 
(4.40) 
(4.41) 



where e is the completely antisymmetric unit tensor of the third rank, (e • n) — e^kUk- 
It is easy to verify that these friction tensors satisfy the symmetry relations [8] 



wrr(l) 
./3a 1 



e?V). (Cm) = £SV> 



a,/5 = l,2,...,iV, 



(4.42) 
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where (j^a\ u )) me ans the transposition of the matrix £pa\w) with respect to the space 
indices i,j = x, y, z. 

Furthermore, the tensors ^p\ijS) and Ca^i 00 ) are independently symmetric both in the 
space indices = x, y, z) and in the indices of spheres a, (3, i.e., 

= Cm- («3%o) t = «3%o. 
G (1) m = e TM, (c (1) h) t = G (1, m- 

By virtue of relations (4.35), (4.36), (4.39) and (4.40), it follows that the components 
of the velocities U p{oj) and £lp{oj) parallel or antiparallel to the vector R a p connecting the 
centers of spheres a and (5 give no contribution, respectively, to the torque and force exerted 
by the fluid on sphere a. 

V. THE SECOND ITERATION 

For the second iteration, the system of equations (2.67), (2.77) is reduced to the form 

E K!iZ» ■ /i 2 L» = - E E • /gL», (s- 1 ) 

Yi 2) H= E e m ■ f*lm{u) = 4tiV3<£ E P? ind (Ra,u;), a = 1, 2, . . . , TV, (5.2) 

m=— 1 /3^a 



where pf >md (R a) uo) is defined by Eq. (2.79) with / fti2m2 (u) — > /^] 2m2 (<^) and the harmonics 

f3,l 2 m 2 I 



/i 1 ) m (a;) are determined by relations (4.4)-(4.12) obtained for the first iteration. Solving 



Eqs. (5.1) and (5.2) up to the terms of the order 6 3 and a 2 [for faf^^i^)] and a 3 [for 
/a^ m2 (^)] and using Eqs. (3.13) and (4.27), we obtain 



= S lt o /S'H + <V + \ 2 f<M( u) + ^ C = * , r, (5.3) 

/2S(«) = m = ± 1. (5.4) 



3 

Here, 
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f(C,2) /. , 5 
f (*.2) /, 

Ja,2ml W J 



f (C ' 2) (, , I 



E ^p( b l) T M(Ra-y,Oj) ■ £ ^(&/ 3 )T M ( J R 7/3 ,Cu) • C/^Cu) + O^ 3 ), (5.5) 

-faP(&a) E T{b 1 )T M (R ai ,uj) ■ E Hip a Pi a iP 9(^)^(^/7/3) ( n 7/3 x ( 5 - 6 ) 



7^a 



£¥7 



-|&{46g5H + («H) T }, 



( = t,r, m = 0,±l, 



6 2 f 



Km' 



m = 0,±1,±2, 



b 2 £ 



K m • £ (l + b^T M (R aj ,Uj) ■ E a /3^0-/3 7 0-7/3 ^(l/7/3) ( n 7/3 X ty^ 
77^a /3^7 



L2t 3 

°aSa \- TV 

/ j 1 ' rami 



15V21 



mi=— 3 



m = 0,±1,±2, 

4cgSH + («H) T }, C = t,r, 
m = 0,±1,±2,±3, 



P (t,i), 



E ( n °/3 ' /ftOO^' 



&S 0) (-) + « 1) M + ^ 2) M, 

7^a /3^7 



b (t,2,l) 



a; 



( = t,r, m = 0, ±1, 
m = 0,±1, 

E • E - i)^ M (^, • u fi (u>) + &^ 



(5.7) 

(5.8) 

(5.9) 

(5.10) 
(5.11) 

(5.12) 
(5.13) 



7^a 



0#7 



>(r,2,l). 



x (l - \K) m E ^Soo • {^V) + (&g?M) T }, - = 0,±1, (5.14) 

= E T a,lmM ' E a P^P (T Pl (T ll3 i>{Vlp) ( n 7/3 X > m = 0, ±1, 

• E( T ( & 7)^(^) - ljaptpafaCTtf ^(2/7/3) ( n 7/3 X + & c 



(5.15) 



E ^aOm( w ) 



7^a 



/^7 



1.2 

15 



x (l - lb«) m E ^So • {4bS 1 ° ) H + K r i 0) (^)) T }, m = 0,±1, (5.16) 



&(C,2,2) ( 

^i,m (' 

l(r,2) 
l a,m 



on 



a; 



!.2 4^ \ - ^i mi 



E ^3 



"a 1 r / / ' ' lm,00 
10 mi=-l 



w) + (fiH) T }» C = *,r, m = 0,±l, (5.17) 



= - E ^ilM • E ( & 7 + h)^T M (R^, to) ■ U p (u), m = 0, ±1, (5.18) 

77^q /3^7 

= E & 7 T Mm( W ) " E a /3^/37 CT 7/3^(l/7/3)( n 7/3 X ^M). 771 = 0, ±1, (5.19) 

7^0 /3^7 

=&s o) M+e ) (-) 
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• 2 (l + ^ + 6 7 )^T M ( J R 7/3 , • C/^H, m = 0, ±1, (5.20) 

= E (! + 6 7) T7 i(^) • E "fi&th^Mirhti x ™ = 0, ±1- (5.21) 

Note that in these relations only the terms up to the order b 3 might be retained. 

Next, we substitute relations (5.5) and (5.6) into Eq. (2.103) and relation (5.7) into 
Eq. (2.104). As a result, we obtain the following expressions for the force and torque exerted 
by the fluid on sphere a for the second iteration: 





N 
0=1 




(5.22) 


F^\u) = 


N 


■ 


(5.23) 


TM(u) = 


N 

-EC?V) 

0=1 


■Uf,(u), 


(5.24) 


T^\u) = 


N 

-ECj 2) H 




(5.25) 



where (cu), A 4 , f = t,r, a,j3 = 1,2, . . . , N, are the friction tensors for the second 
iteration, which can be written in the form 

£f(u) = UpP(b a )p(bp) E p(b^%T M (R ai , u) ■ T M (Rfr, u), (5.26) 

= -^papp{K)q{b p ) r(b 1 )a M a lP i)(y M )T M (R ai ,uj) ■ (e • n Pl ) , (5.27) 

£$ 2) M = £a£/3a a p(bp)q(b a ) ^ r^a^a^^y^) (e • n a7 ) • T M {Rfr, w), (5.28) 

7^a,/3 

3 (X 2 

C/3 2) (^) = t(— ) C a q(ba)q(bp) ]T xfa)<^<Tfratf$(ycnMy(h){(ncn • nfh) 1 

- n^n^}, (5.29) 

where 

xW = l + &+^ 2 + y. (5.30) 
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Using the explicit form (4.24) for the modified Oseen tensor T M (R a /3,u), it is useful to 
represent the tensors T M (R ai ,uj) ■ (e ■ n^) and (e ■ n aj ) ■ T M (Rp 1 ,uj) in relations (5.27) 
and (5.28) as 



T M (R cer ,u) ■ (e • n Pl ) = 



(e • n Q7 ) • T M (R /3l ,uj) = 




In the particular case of two spheres (a,/3 = 1,2), according to relations (5.26)-(5.29), 
all friction tensors $j£p 2 \w) = 0, [i, v — t,r, for (3 ^ a, i.e., these components are caused 
only by three-particle interactions. 

Just as the friction tensors ^^{oj) for the first iteration, the friction tensors ££p 2 \w) 
for the second iteration satisfy the symmetry relations [8] 

(CV)f = ev), (c (2) m) t = a (2) H, (c } m) t = cvx 

a,(3 = 1,2,... ,N, (5.33) 

Substituting relations (2.99), (3.40), (4.34), and (5.22) into (2.98), relations (3.41), (4.35), 
and (5.23) into (2.100), relations (3.46), (4.36), and (5.24) into (2.101), and relations (3.47), 
(4.37), and (5.25) into (2.102), we obtain the following expressions for the forces and torques 
exerted by the fluid on the spheres up to the terms of the order b 3 : 





N 
(3=1 


■Uf,(u), 


(5.34) 




0=1 


■ 


(5.35) 


t2>(") = 


13=1 




(5.36) 


T t\u) = 


N 






13=1 




(5.37) 
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Here, 

CM = Uu)8*P 1 + ^A^H, a, (3 — 1,2, . . . , N, (5.38) 

CM = &{u>)5af, I + a, (3 — 1,2, ... ,N, (5.39) 

ip{y ai )T M {Raj, v) • (e • n Q7 ) , (5.40) 

OH = -^a^g(6^)|r(6 a )(T a/3 (T /3a ^(t/ Q/ 3) (e • n a/3 ) 

+ £aP{b a ) T(6 7 )^ 7 f7 7 ^V(Z//3 T ) T M(-Ra7, ^) * ( e ' n /3 7 ) \ , P (5-41) 

Ca(^) = £la a p(b a )q(b a ) ^(b^a^a^iy^) (e • n a7 ) • T PI (R a7 , u), (5.42) 

CjsM = -^a a g(6 a )|r(6 /3 )o- a/3 or /3a V ; (ya/3) (e • n Qj9 ) 

- £/3P(M XI T (b 1 )(J ai a ia 'ip(y ai ) (e • n a7 ) • T M {Rp 1 , (3^ a, (5.43) 

*a» = £aP 2 (U J2p(b,)^(T M (R ay ,iu)J, (5.44) 



7^a 



AjsM = -Zppib^pib^TMiR^u) 

- J2 P(bj)^T M (R ay ,uj)-TM(R^,uj)\, (3^ a, (5.45) 
3 

^ r Q (w) = ^ q 2 (b a ) x{b 1 )(rl 1 a lol i) 2 (y ai )(I - n Q7 n a7 ), (5.46) 

7^a 

3 o 2 



A^(w) = q{b a )q{bp)\ yl p \h 2 (y aP )(I - 3n afS n a/3 ) + 2h (y afS )I^ + 



4 \a a J 



x S x( & 7) CT a 7 CT /37 CT 7/3V ; (l/a7)V ; (l//3 7 ){(n a7 -n / 3 7 )/-n /37 n Q7 U, /3 ^ a. (5.47) 

Note that in relations (5.41), (5.43), (5.45), and (5.47), the terms proportional to J2 
correspond to the three-particle interactions and, hence, they are absent in the particular 
case of two spheres (a, (3 = 1,2). 

Note that the friction tensors for the first and second iterations satisfy the symmetry 
relations given by Eqs. (4.42) and (5.33), and, hence, these relations also hold for the total 
friction tensors £^(w), v — t,r, i.e., 

a,P= 1,2,..., N. (5.48) 
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As was mentioned above, all required quantities (including the friction tensors) are ob- 
tained without imposing any additional restrictions on the dimensional parameter y a p. In 
particular cases, the expressions for the friction tensors can be essentially simplified. We 
consider the following two opposite modes: Rey a/ 3 ^> 1 and \y a p\ <C 1, which is typical of 
the dilute suspensions and colloidal crystals, respectively. 

For the first mode, the frequency belongs to the range 

v 2u 

> M > 7^, (5-49) 



"'max v mm 



where a max is the radius of the largest sphere and R m in is the distance between the centers 
of two closest spheres. 

Then we can use the asymptotics 

h n {z) ~ exp(-z), n = 0,1,2,..., (5.50) 
Zz 

valid for large z, which enables us to simplify the friction tensors (5.38)-(5.43) to the form 
£%3( u ) ~ £a(u)5 a/3 I + ^ £aP(b a )p(bp)\ (1 - <M (i - 3n a/3 n a/3 ) 

1 I Vaf3 

+ f P(b-y) -Y—^ {l + 9(n Q7 • n f3l )n ai n l 3 1 - 3n ai n ay - Sri^n^} 1, 

a,P=l,2,...,N, (5.51) 





a, (3 = 


1,2,. 


..,7V, 


(5.52) 


CM « o, 


a, /3 = 


1,2,.. 


,.,7V, 


(5.53) 


ejH « o, 


a, f3 = 


1,2,.. 


.,N. 


(5.54) 



This mode is characterized by an exponentially small coupling between the translational 
and rotational motions of the spheres proportional to exp(— y a g) for f3 ^ a as well as between 
the rotational motions of different spheres. Taking into account that in the stationary case, 
O' O' and for ^ ^ a ' are P ro P ort ional to V^a/j and 1/^/3, respectively, [29], 
we may conclude that the retardation effects lead to an essential decrease in hydrodynamic 
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interactions between the spheres. For the mutual translational friction tensors, this influence 
is not so considerable because the behavior of as 1/R a p in the stationary case [4] is 
replaced by 1/R^p in the nonstationary case for large distances between the spheres. 
The second mode, \y a p\ <^ 1, corresponds to the low frequency domain. In this case, 

^(yap)^l- y -f+ V -f + 0(yt,), (5.55) 

T M (R a p, u) » T M {R aP ) - — -^g — / + 0(y 2 a/3 ), (5.56) 

where 

T M (R a p) = — (j + n Q/3 n Q/3 ) + (a 2 ap + <r| Q ) Q I - n a/3 n a/3 ) | (5.57) 

is the modified static Oseen tensor [42]. Note that within the framework of the second 
iteration, only terms proportional to 1/R a p should be retained in (5.57), i.e., 

T M (R a p) ~ \ (i + n a/3 n a(3 y (5.58) 

Substituting (5.55) and (5.56) into (5.38)-(5.47) and retaining only the main frequency- 
dependent terms, we represent the friction tensors in the form 

^H = ^ + ^H. H," = t,r, a, (3 = 1,2,..., TV. (5.59) 

Here, fi,v = t,r, a, (3 — 1, 2, . . . , N, are the stationary friction tensors determined up 
to the terms corresponding the second iteration [4,27,30,39]: 

£p = (Sapl + Kp) , a, (3 = 1,2,..., N, (5.60) 

Q, = £ (<W J + Kp) , a, (3 = 1,2,..., N, (5.61) 
3 

iaa = -i a a a - vl^ja (e • n a7 ) , (5.62) 

Co = — ^aa) (5.63) 

= Spap^apvpa (e • n a/3 ) 

+ 7 H o^ipopi {(e • n^) - n Q7 (n a7 x n^)} 1, /3 7^ a, (5.64) 
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Caf3 — —^aaa^CapVpa (e • Tl a p) 

~ 7 Z o-^a^a^^e ■ n Q7 ) + {n M x n ai ) n^j 1, (3 ^ a, (5.65) 
9 

X aa = E ( J + 3n Q7 7l Q7 ) , (5.66) 
X % = -^|^/3a (i" + ™ Q/3 n Q/3 ) 

-| E (7^/3^(1 + n ai n ai ) ■ (I + np 7 np 7 )\, P ^ a, (5.67) 
3 

^«« = |E ^7^" ( J - ^7 n «7) , ( 5 - 68 ) 

Kp = i 1 ~ 3n a pn a p) 

+ | (— ) Z CT /37 a 7/3 a a 7 {( n «7 •n /37 )/-n / 3 7 n a7 } 1, /3^<x (5.69) 



2 Va a / 7 ^ Q ,/3 

The frequency-dependent terms A^(u), /i,i> — t,r, a, (5 — 1, 2, . . . , N, which vanish as 
uj — > 0, have the form 

{3 3 f 

6gl - (l - <W) - ^(o- Q/3 + op^j (i + n a/3 n a/3 ) - - ^ o- /37 6 7 (l 

+ npjUp^ + <J ia bp(l + n Q7 n Q7 ) |, (5.70) 

O) « £ { f (<W - I + (l - <W) f ^ (J + n a/3 n a/3 ) + | (£) 2 

X S {2^7^/37^7/3^7 - (yfoVlftf* + OA) ( X + 6 7) + -( ^^ 6^ 

7^q,/3 k ' 6 °"7 

+ ^I6j6 7 ) J|(n Q7 • n^ 7 )l - n^n^J j, (5.71) 

^JsM ~ ^a^aj (l - <W) (y - o" a /30-/3a) (e • n Qj g) + ^ 073 7 o- 7/3 (e • n^) |, (5.72) 

Aapi^) ~ CaQa^l (l - <W) (-^ - VapapJ) (e • n a pj - a ai a ia (e ■ n a7 ) 1. (5.73) 

According to (5.70)-(5.73), in the low-frequency range, the main frequency-dependent 
components of the friction tensors are proportional to uo 1 ! 2 . For the rotational and coupled 
friction tensors, we also retain terms proportional to uj (for the rotational friction tensors, 



41 



even terms proportional to uj 3 / 2 ) because they are of the order a n [n = for A r ^ a (uo) and 
A^Ju) and A r *Ju) for (3 ^ a and n — 1 for A r ^{u) for /3 ^ «], whereas the terms 
proportional to uo 1 / 2 are of the order of a m [m = 2 for A^^uj) and a, (5 — 1, 2, . . . , N, 

and m = 4 for for a, /3 = 1, 2, . . . , N] with m > n. 



VI. MOBILITY TENSORS 

Substituting Eqs. (5.34) and (5.35) into (2.96) and Eqs. (5.36) and (5.37) into (2.97) 
and solving the resulting system of equations for the quantities U a (cu) and Q a (cu), we can 
represent these quantities as follows: 

N 

u*{u) = = -£ {aOM • + • , (6-i) 

0=1 
JV 

n«H = = - E K» • + ■ T p (u,)} , (6.2) 

a=\ 

where /j^Jcu) and ^JuS) are, respectively, the translational and rotational mobility tensors, 
and n^piui) and //^(u;) are the mobility tensors that couple translational and rotational 
motions of the spheres. Up to the third order in the parameter b and the same orders in 
the parameter a that for the corresponding friction tensors, i.e., up to terms of a 2 , a 4 , and 
a 3 , for the translational, rotational, and coupled mobility tensors, respectively, the final 
expressions for the mobilities can be represented in the form 



= tfMl + J- E b 2 a ay a y Jexp(-2y a ^)I + - (j + 3n a7 n Q7 ) 

exp(-y Q7 )(j - 3n Q7 n Q7 )J, (6.3) 



2 

2 



A*$jM = {l + - [&* (1 - 6 a ) + $(1 - M] }r M (^, + — ]T & 7 ^ 7 ^|exp(-y a7 ) I 

- (I - 3n a7 n a7 ) } • {exp(- y/37 )J - (/ - Sn^n^) }, ^ a, (6.4) 

Maa(w) = + 777 E & ? "a 7 CT 7"^ 2 (^7) ( J ~ ™a 7 ™a 7 ) , (6.5) 

A*I£M = 3^r{^2(ya/3)(3n Q ^n a/3 - /) - 2/i (z/a/3)J r } 
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+ 7r^(— ) 2 12 b 2 y a ai a ia a^(y ay )^(y^){(n ai -n^I -n^n ai }V ^ a, (6.6) 
1 



(6.7) 



6 



x ||exp(-y QT ) - w ^" 7 ^ |( e • n^) - ^ w{y ai )n ai (n ai x n^) j, /3 7^ a, (6.* 



Maa(w) = (Maa(^)) T = — E ^^F^(l/a7){ eX P(-^7) " 

O071V/ ^ q /t a7 L 

'a/3) 



w{y, 



'07 ; 



e ■ n 



cry 1 5 



(6.9) 



X 



|exp(-y^ 7 ) - w(kV ^ |( e • n Q7 ) + jj w{y Pl )(n Pl x n Q7 )n^ 7 |, /3 ^ a. (6.10) 



Here, 



A4M 



Co 



Set 



6 2 

l- f (l- 6 . 



(6.11) 
(6.12) 



are, respectively, the frequency- dependent translational and rotational mobility tensors of a 
single sphere calculated up to the terms of order b 3 [31], and 



6 4 ~ 
w(x) = — xh 2 (x), 



(6.13) 
(6.14) 



Note that the obtained mobility tensors defined by relations (6.3)-(6.10) satisfy the 
symmetry relations [8] 

(m£,M) t = /*$,M, (/*&M) r = /W> «M) T = ^H, «,/5 = i,2,...,iv. 

(6.15) 

In [31], the frequency-dependent mobility tensors are derived up to the terms of the third 
order in two parameters b and a, i.e. up to the terms of b n <7 m , where n + m = 3 ("n + m = 3 
approximation"). According to one of the main conclusions in [31], the frequency-dependent 



43 



mobility tensors determined in this approximation contain only the terms corresponding to 
the two-particle interactions, which is similar to the behavior of the frequency-independent 
mobility tensors determined up to the terms a 3 ( or <r 4 for fJ^p) [4,29]. Moreover, the diagonal 
(f3 = a) translational and rotational mobility tensors determined in [31] are independent of 
the hydrodynamic interactions between the spheres, i.e., 

M*» = /C» = (6.16) 

and, hence, the self-interaction of a sphere due the action of the fluid induced by this sphere 
and reflected from the rest spheres is not described within the framework of approximation 
n + m = 3. Furthermore, there is no coupling between translational and rotational motions 
of the same sphere, i.e., 

/O)=0, /&M = 0. (6.17) 

Unlike [31], expressions (6.3)-(6.10) represented in the form of expansions in the param- 
eters b and a contain the terms up to b 3 a m , where m = 2 for //^(u;), m = 3 for /^(u;) 
and Hapiou), and m = 4 for fjt^(u), i.e., n + m > 3. An account of such terms in the 
mobility tensors leads to some new effects, caused by the contribution of the three-particle 
hydrodynamic interactions to the mobility tensors. This contribution is described by the 
terms proportional to the sum in Eqs. (6.4), (6.6), (6.8), and (6.10). It follows that the 
main terms corresponding to the three-particle hydrodynamic interactions are proportional 
to b 2 a 2 for fj^p(u), b 2 a 4 for fj£p(u), and b 2 a/R 2 a p for fj^p(u) and fj, r Jp(u). In addition, the 
terms proportional to i n the diagonal mobility tensors ^ a {oo) and fi r J a (uj) defined by 
relations (6.3) and (6.5) describe the self-interaction of sphere a due to the action of the 
fluid induced by this sphere and reflected from the rest ones. This effect is of the order of 
b 2 a 2 for ^ a {oS) and b 2 a A fj, r J a (uj). According to Eqs. (6.7) and (6.9) for the mobility tensors 
^ r a (w) and /J>a a (u}), there is a coupling between translational and rotational motions of the 
same sphere, which is of the order of b 2 a/B? . 

Now, let us examine the mobility tensors defined by Eqs. (6.3)-(6.10) in the approxima- 
tion n + m = 3. We see that the relations for /i" a (u;), ^ a {uj) and fj,^p(u), n r *p{uj), where 
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a, (3 — 1, 2, . . . , N, coincide with the corresponding relations given in [31]. For the rotational 
mobility tensors fjt^^u), where (3 ^ a, we get 

A*E/>M = ~2^j{^ y<x ^( I ~ 3na/3na/3 ) + 2 ^o(ya/3)i"|, (6.18) 

which differs from relation (4.5) in [31] by the second term in the braces. [In terms of the 
notation used in the present paper, relation (4.5) in [31] is reduced to the form (6.18) with 
the substitution of (l/2)h (y a p) (jl - 9n Q/3 n a/3 ) for 2h (y al3 )I.] 

The translational mobility tensors fi^^u), where (3 ^ a, in approximation n + m = 3 is 
reduced to the form 

^(u)=T M (R a p,u) + l(b 2 a + b 2 p )T M (R a p), (6.19) 

which essentially differs from the translational mobility tensor given by Eq. (4.3) in [31] 
because of the presence of the terms necessary for holding the symmetry relations [8]. 

In the particular case of two spheres (a, (3 = 1,2), Eqs. (6.4), (6.6), (6.8), and (6.10) for 
(3 7^ a are simplified to the form 

M2M = ft"H = {l + I [&i(l - 6i) + b 2 2 (l - b 2 )]^T M (R 12 ,u;), (6.20) 
= = -^fy { h ^){' - 3 "i2"i2) + 2Ao(i/i 2 )/}, (6.21) 

*> = S;( 1+ t)« ( e -4 (6 - 22) 

and the mobility tensors /^(u;) an d M^U^) are defined, respectively, by relations (6.22) and 
(6.23) by interchanging the indices 1 and 2. 

Thus, in the case of two spheres, the translational and mobility tensors are independently 
symmetric in the space indices = x,y,z) and in the indices of spheres a, (3 — 1,2, i.e., 

= A*"a( w )» (O)f = mSjH, a,/3 = 1,2, 
/x^M = /^», (k») T = /x^H, a, 0=1, 2. 
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Analogously to the friction tensors, we investigate the behavior of the mobility tensors 
in two special opposite cases. First, we consider the case Rey a p 1 corresponding to large 
distances between the spheres. Using asymptotics (5.50) and 

6 

w(Vafi) ~ — , ( 6 - 25 ) 

we reduce the mobility tensors (6.3)-(6.10) to the form 
li%{u>) » fia^Sap <W) ^ {l + | [ft* (1 - 6 Q ) + b}(l - bp)] } (3n a/3 n a/3 

(7 (X 3 1 

- I) + XI " 7 2 7/3 {/ + 9(n Q7 ■ npjjnajnpj - 3n a7 n Q7 - Sn^n^j L 

a, /3= 1,2,..., TV, (6.26) 

aC/jM ~ A*»<W J, a, = 1, 2, . . . , N, (6.27) 

A*^H»0, a, /3= 1,2,..., TV, (6.28) 

A*3jH»0, a, /3= 1,2,..., TV. (6.29) 

Note that the nonstationary mobility tensors fj^^u) and a, j3 = 1, 2, . . . , N, and 

the mutual rotational mobility tensors P ^ a, exponentially vanish as exp(— y a p)i 

whereas in the stationary case, for j3 ^ a, these quantities decrease as l/R^p an d l/R^p, 
respectively, [29]. 

Relations (6.27)-(6.29) coincide with the corresponding relations given in [31], whereas 
the translational mobility tensor defined by (6.26) differs from the corresponding result given 
in [31]. Indeed, in approximation n + m = 3, relation (6.26) is simplified to the form 

l r 2 i 
»%{u)^va{u)8 a pl- —— l + -(b 2 a + b£) (i-Sn^nJ), (6.30) 



which differs from Eq. (4.6) in [31]. Indeed, in terms of notation used in the present paper, 
the latter takes the form 

1 A , 2, 2 aa 



A*$jM ~ A*aH<W 1 ~ 7 p 2~ ( 1 + q ) ( J ~ 3n «/3 n «/3) > ( 6 - 31 ) 
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which does not satisfy the symmetry relation [8]. 

Finally, consider the mobility tensors in the low- frequency domain \y a /3\ <C 1. Using 
Eqs. (5.55), (5.56) and the series 



"(Vctf) « 1 - ^ + ^!& + -, (6-32) 



and retaining only the main frequency- dependent terms, we simplify the set of the mobility 
tensors given by Eqs. (6.3)-(6.10) to the form 

M^H=M^ + ^H, W = t ^ a,P = l,2,...,N. (6.33) 

Here, /u^, //, z/ = t,r, a, (3 = 1,2, . . . , N, are the stationary mobility tensors calculated 

within the framework of the second iteration [4,29,39]: 

If 3 I 

Ma/3 = |^/3 J + (! - <M ^ ff a/3 (I + n a(3 n af3 ) J , (6.34) 

Ma/3 = ^ {<W-T + ( x ~ M (3ri a/3 n a/3 - J) I , (6.35) 

Ma/3 = (1 - M U ( e ' > ( 6 - 36 ) 

Ma'/3 = Ma/3- (6-37) 

The frequency-dependent terms M^(u), \i,v = t,r, a, (5 — 1, 2, . . . , N, which vanish 
as uj — > 0, have the form 

Ml»«-^J, (6.38) 
6 2 3 

+ ^r(— ) 51 b^a ai (x ja a 2 ^Un aj ■ npjjl - n Pl n ai \ (6.39) 



+ E ^ P^P |( e ' ^ - ™a 7 (™«7 X n /3 7 ) } f > ( 6 -40) 
7 ^a,/3 - n -a7- n -/3 7 I- I 



El2 °7Q 
7 D p 

7 ^a,/3 - n -a7- rl /37 



- E 6 7 fl^fl {( e ' n «0 + ( n /?7 X ^a 7 ) n /3 7 } |- ( 6 -41) 
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The frequency-dependent component M^^uj) of the translational mobility tensor co- 
incides with the corresponding result given in [31]. In [31], the mobility tensors in the 
low-frequency range are calculated to the first order in the parameters b and a, and, within 
the framework of this approximation, for (3 ^ a, 

A*3,M = 0, A*^H = 0, a, (3 = 1,2,..., TV, 

= 0, (3 ? a. (6.42) 

As is seen from relations (6.39)-(6.41) derived in the present paper within the framework 
of the approximation discussed above, the mobility tensors //^(u;), /x^g(cj), and //^(a;) for 
a, (3 — 1, 2, . . . , N, are not equal to zero but proportional to u in the low-frequency range. 

VII. CONCLUSION 

In the present paper, we calculate the distributions of the velocity and pressure fields in 
an unbounded incompressible viscous fluid induced in the nonstationary case by an arbitrary 
number of spheres moving and rotating in it as well as the time-dependent forces and torques 
exerted by the fluid on the spheres. All the quantities of interest are expressed in terms of 
the induced surface force densities without imposing any additional restrictions on the size 
of spheres, distances between them, and frequency range. Within the framework of this 
approach, we derived the general relations for the translational, rotational, and coupled 
friction and mobility tensors in the case of an arbitrary number of spheres up to the terms 
of the order b 3 and a 2 , a 4 , and a 3 , respectively, and analyzed the contributions to the 
mobility tensors caused by the effects of the three-particle hydrodynamic interactions and 
the self-interaction of spheres. We investigated as well the influence of retardation effects 
on the behavior of the friction and mobility tensors in the low-frequency range and for large 
distances between the spheres. 

Finally, we formulated the procedure which allow one to calculate both the velocity and 
pressure fields in the fluid and the forces and torques exerted by the fluid on the spherical 
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particles up to the terms of the order a p , where p is any positive integer. With the use of 
this procedure, the hydrodynamic interaction in a viscous fluid can be described analytically 
within to the terms of any given order in the parameter a. 



APPENDIX A 



The quantities K l ™ M and W^ ulm , where h, l 2 , 1 > 0, h > m 1 > ~h, l 2 > m 2 > -l 2 , 
and I > m > —I defined by relations (2.50) and (2.54), respectively, can be represented 
in the explicit form in terms of the Wigner 3j-symbols [41,44] as follows (for details, see 
[39,42]): 

7 

-^(2h + l)(2/ 2 + 1)(2Z + 1) (-l)™i+™ < 



jh-h+l 

Ts-hm2 _ _ 

l\m\,lm 3^/7T 



where 



-( 


-i; 


m 1 
2V 


Uj- 




(2 


I 


L \ 


X 










\ k 


—m 


m — k 



3 — Jmin 



2 / L 




l\ l 2 I 




h h L 




—mi m,2 m 



( 



\ 1 



l\ l 2 L 
^ —m\ m 2 k — m j 



h,k,l> 0, 



L — 1 + 2j, 



1, if / = 
0, if 1 = 1 
-1, if I > 2, 



(Al) 



and 



K 


= V 3 ^ exex ~ e y e y 2 e z e z) = \i 


/| (eie_i + e-id + 2e e ) , 


K 1 




-V2 (e_ie + e e_i) , 


K 


= K] = - V2 (eie + e ei) , 




K, 




2e_ie_i, 


K 2 


= K* 2 = 2e iei , 





(A2) 
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l\m\,lm 



2^ 



(2Z 1 + 1)(2Z 2 + 1)(2/+1)(-1)"™ E 4 E (2^ + 1) 



fc= — 1 J=ir7 



( I I L \ ( U U l \ ( I I L \ ( h l 2 L 



x 







'1 ^2 





fc — m m — k 



-mi m 2 k — m 



J 



h,l 2 ,l> o. 



(A3) 



where 



L = l + 2j-l, 



Note that K\\Z\ M * only for 



1, if I = 
0, if Z > 1. 



Z = Zi + Z 2 - 2p > 0, 



(A4) 



where p = -1,0, 1, . . .,p m ax, Pmax = min([(Zi + Z 2 )/2] , 1 + min(Zi, Z 2 )), H is the integer part 
of a, and min(a, 6) means the smallest quantity of a and 6, and Jm ^ only for 



/ = h + l 2 - 2p + 1 > 0, 



(A5) 



where p = 0, 1, . . . Pmax = min([(Zi + Z 2 + l)/2] , 1 + min(Zi,Z 2 )). 

This means that the quantities i 7 } 1 ; 2 ,/(r a , ag, i? Q( g, a;) and C; 1 « 2i ;(r Q ,, ag, i2 aj g) defined by 
relations (2.49) and (2.53) for (3 ^ a should be determined only for the values of / given by 
(A4) and (A5), respectively. The explicit form for these quantities is as follows [39]: 

^3/2 r (h + i 2 + 



2k 



F h i 2il (r a ,a p ,R af3 ,uj) = (-l) p — lj h {x a )ji 2 {b p )hi(y a p)-5 hh+ i 2+2 , x , 

n V \ 2 r(Zi + |)r (l 2 + | J 



y a l3 n al3 



Chh,l( r ay a /3, Rap) — ^Ml+k+l 



x /- h + h-2p>0. p= -1.0.1 



(A6) 



r(z 1 + z 2 + |) ^ 
2 r(z 1 + |)r(z 2 + l) i^+ 2 ' 

Z = Z 1 + Z 2 -2p+l >0, p = 0,l,...,p maa; , (A7) 



for r a < R a /3 — and 
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2K [~ ~ ~ 7T 3/2 

F hh j(r a ,ap,R a p,uj) = (-l) i2 ~ p — \ h h (x a )j h (b )ji(y a p) - <5 Ml _ i2 _ 2 — - 

711] [ Z 



Til 1 1 



x 



1 ' 9. 



h r>h-2 



r(z 2 + |)r(z 1 -z 2 -i) J' 

l = h+l 2 -2 P >0, p= -1,0,1,..., p m a X , (A8) 



^hhA r oc, 0>/3, Map) - —— 7 a \ „/ 1+1 ^ ^,/i-/ 2 -l 



2 r(/ 2 + |) rfc 



r(z 1 -z 2 + |) 



+ 5« 



'2— 3 



X 



r(z 1 -z 2 -|) 



(7 



/3a 



1 



2Zi -1 V^afl/ 2/ 2 + 3 2/ + 3 



/ = Z 1 + Z 2 -2p+l > 0, P = 0,1, ■■■,Pmax, (A9) 



for r a > R a/3 + a/3. 
Here, ji(x) = 



<n/(2x)I l+ i(x) and = \Jtt / (2x)K l+ i(x) are the modified spherical 

Bessel functions of the first and third kind, respectively, [43], T(z) is the gamma function, 
<rpa = ap/R a p, %a = Kr a , and y a p = nR aP . 

In the particular case Z = 0, the general relations (Al) and (A3) for the quantities 
K\IZm and W\\Z\ M are simplified to the form 



ih—h 



and 



'2 /■ 

= ycZfifc (J - n fe n fc ) Yi imi (6 k , ip k )Yi 2m2 (9 k , (p k ) 



3^ 



<W 2 * 



mi,m2 



/ 3(2Z 1 + 1)(2Z 2 + 1) ^ 2( ^ 



mi 



/ /i Z 2 2 ^ 







X 6 k 
k=-2 



W 1 ' 2 " 12 - - 





( k 


Z 2 2^ 


> 


mi 




m 2 k j 


l< > 




-mi 





(A10) 



jh— li— 1 



2v^F 



=- y*rffi fe n fc ^* mi (0 fc ,^ fc )y j2ma (0 fc ,^ fc ) 

/ /i z 2 1 ^ 




^(2Z 1 + l)(2/ 2 + l)(-l) mi 



V 



i 

X ^ ' 3k,mi—m2 
k=-l 



-mi m 2 k 



(AH) 
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With regard for the properties of 3j-symbols, we have 

0,1, if/i = 0,l 



KS, oo ^0 only for Z 2 = h + 2p > 0, where p 
and 



(A12) 
0,±1, if/i>2, 



, 0, if h = 

Wft oo ^ only for Z 2 = Zi + 2p + 1 > 0, where p = \ (A13) 

0,-1, if h > 1. 

According to (A13), the quantity Pj 1 2,j(r Q , a^, i? Qi g, a;) contained in relation (2.40) with 
Z 2 = 2 should be determined only for Z = Zi ± 1 > 0. For this value of Z, we have [39] 

2 ~ 

^i2,/( r a,a/3,-Ra/3,^) = (-l) p — 3h{ x a )j2{bp)h l (y a(3 ) (A14) 
for r a < i? Q/ 3 — ag and 

2 ~ 

Ph2,i( r oc, a i3, R ai3,u) = (-l) p+1 — h h {x a )j 2 {bp)ji{y a p) (A15) 

7T77 

for r a > it! a/3 + ag. 

According to (A12) and (A13), the quantities F^; 2 (r a , a a , uj) and Ci^ij-a, a a ), should be 
determined only for Z 2 = h + 2p > 0, where p = 0, ±1, and Z 2 = Zi ± 1 > 0, respectively. As 
a result, for r a > a a , we get [39] 

F hh (r a ,a a ,uj) = (-l) p — lh h (x a )j h (b a ) - 5 l2th _ 2 (k ~ 7j) -3 (— ) } » 

Z 2 = Z x + 2p > 0, p = 0,±l, (A16) 

C h i 2 (r a ,a a ) = 5, 2)Jl _i— ( — ) + <5 ra , a — (<J{ 2) i 1+ i - <^i-i) , h = h ± 1 > 0. (A17) 
For P 12 (r a , a a , w), where r a > a a , we have 

2 ~ 

^1,2(^0, a a ,w) = — h 1 (x a )j 2 (b a ). (A18) 

7T?7 

In the particular case r Q = 0, we have 

a' 2 

C llh , l (0,ap,R a p)=6 ll>0 6 l>h+1 -£-, Z = Z 2 ±1>0, (A19) 

lim C hh (r Q , a a ) = <J Jl>0 <^ 2 ,i 4" ■ ( A20 ) 
ra—o at 
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APPENDIX B: ON THE VALIDITY OF SOLUTIONS OF THE LINEARIZED 

NAVIER STOKES EQUATION 



All relations obtained in the present paper are true within the framework of the linearized 
Navier-Stokes equation. Using the explicit form for the velocity and pressure fields of 
the fluid induced by a single sphere or a system of noninteracting spheres, we consider 
several problems of the correctness of the linearization of the Navier-Stokes equation in the 
nonstationary case. The Fourier transform of the linearized Navier-Stokes equation (2.1) 
has the form 

Vp(r, u) — {iup + r)A)v(r, u) = — pVy?(r, u), (Bl) 

where A is the Laplace operator. The linearized equation (Bl) is valid [8,11-15] provided 
that the force acting on a unit volume of the fluid due to the fluid viscosity 

F vis (r,t) =r]Av(r,t) (B2) 

is much greater than the nonlinear (inertial) term 

F nl (r,t) = -p(v(r,t)-V)v(r,t). (B3) 

To estimate the possibility to neglect the nonlinear term, we use the following condition: 

\F vis (r,t)\ > \F nl (r,t)\, (B4) 

where the bar in (B4) means averaging over a certain time interval. 

Since in the approximation of noninteracting spheres, the rotation of spheres does not 
make a contribution to the fluid pressure [relation (3.21)], we rewrite Eq. (Bl) in terms of 
the notation used in the present paper as follows: 

Vp mind (r, u) - {loop + 77A) (v mind (r, u) + v^ r ' 0)ind (r, u)) = 0, (B5) 

where 
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N 

v^ ind (r, u) = £ vf mind (r, w), (B6) 

(3=1 

N 

v (r,0)ind( r ^ ^ = v {S,r,0)ind( r ^ ^ + w (V)ind( r> ^ = £ ^MW^ ^ ( B7 ) 

0=1 

p(t,o)<nd (r> w) = £ p^°) <nd ( r> w ). (B8) 
/?=1 

We represent the Fourier transform of the viscous force F ms {r,uj) as a sum of two 
components 

F ms (r,cu) = F vis ^(r,cu) + F vis{r \r,cu), (B9) 

where 

F vis{t \r,u) = r]Av^ md (r,cu), (BIO) 
F vis(r) {r, u) = r]/\ v ^ ind (r, u) = -iupv^ md {r, u) (Bll) 

are the viscous forces due to the translational motion of the spheres and their rotation, 
respectively. As is seen from relation (Bll), the nonstationarity of the problem leads to the 
appearance of the component of the viscous force induced by rotation of spheres. 

First, we consider the case of a single sphere a. To verify condition (B4), the general 
relations for the fluid velocities v^ ind (r,u) and v^ md {r,u) in the entire frequency range 
are necessary instead of particular relations (3.29) and (3.25). To this end, we use the well- 
known results for the fluid velocity induced by a single sphere moving with arbitrary velocity 
U a (t) [16] 

3 a n 



v mind M = ^^{p(6 a )(3n a n a - /) + {(l + x a + x 2 a )l 

- (3 + 3x a + xty n a n a } exp(-(x a - b a )) }u a (u) (B12) 

or rotating with angular velocity Q a (t) [14] 

(r, w) = (—) \ + Xa exp(-(x a - b a )) (p a (u) x r a ) . (B13) 



v (r,0)ind( 
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Substituting (B12) and (B13), respectively, into (BIO) and (Bll), after certain trans- 
formations, we obtain the following relations for the Fourier transforms of the translational 
and rotational components of the viscous force: 

F^(r^)= V ^§{(l + x a + xl)l 

- (3 + 3x a + 4) n Q n a | exp(-(x a - &«)) U a (u), (B14) 
F ms(r \r,tu) = 77% I±|* exp(-(:r a - & a )) (n a (w) x r a ), (B15) 

which are valid for any r a > a Q . 

It should be noted that the first term in (B12), being the main term of the fluid velocity 
induced by the moving sphere for distances Rex a ^> 1, gives no contribution to the viscous 
force. As a result, the quantity F ms ^\r,u) exponentially decreases in this range despite 
the fact that v^ md {r,u) is proportional to 

Lets us investigate the validity of solutions (B12) and (B13) of the linear Navier-Stokes 
equation in two particular cases corresponding to a sphere oscillating in a fluid with trans- 
lational U a (t) = U a cos out or rotational ft a (t) = ft a cosut velocities, where the frequency 
of oscillations uj is such that \b a \ <C 1. We investigate these problems in the domain far 
from the sphere: Rex a 3> 1. To this end, using (B14) and (B15), we represent F ms ^\r,t), 
where ( = t, r, as follows: 

F**t\r,t) = Be(^ 8 «\r,u)&v(-iut) > ), ( = t,r, (B16) 

where 

F^\r,u) » ^exp(-x Q )(j- n a n a ) U a , (B17) 

F ms ^(r,uj) » r ? S exp (_ Xa ) (n Q ( W ) x r a ). (B18) 

Substituting the fluid velocities defined by relations (3.29) and (3.30) for Rex a 3> 1 into 
(B3) and averaging the estimated forces, we obtain that in this space domain, inequality (B4) 
is true for a sphere rotating with oscillating angular velocity for any distances r a (except 
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for the particular case where r a || O q (c<j) for which both F ms ^ r \r,u) and the corresponding 
nonlinear force are equal to zero), whereas, for a sphere oscillating with translational velocity, 
inequality (B4) is reduced to the form 

/ 3 \ 5/2 a / \l/2 

| sin^| exp(- Rex a ) > R ( — — - ) \p{b a )\ 2 — (sin 2 ^ + 8 cos 2 ^(1 + cos 2 V) J , (B19) 

\2\r a \ z / r a v ' 

where 

R = ^ (B20) 

is the Reynolds number, U a = \U a \, and ip is the angle between the vectors U a and r a . 

Since Rex a ^> 1, inequality (B19) is not valid for any r a in this domain. Therefore, 
expression (3.29) for the fluid velocity induced by an oscillating sphere far from it obtained 
within the framework of the linear approximation cannot be regarded as correct in this space 
domain for the problem considered because, in the general case, the nonlinear term is much 
greater than the viscous force. 

Using the results obtained above, we can make analogous conclusions for a system of 
oscillating spheres in the approximation of noninteracting particles. In particular, the fluid 
velocity and pressure induced in the far zone by spheres oscillating with velocities U pit) = 
U i3 cosut, where (5 = 1, 2, . . . , N, should be determined using the nonlinear Navier-Stokes 
equation because the linear solutions defined by relations (3.33) and (3.34) contradict the 
required condition (B4) of linearization of the Navier-Stokes equation at any point in this 
domain. 
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